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ABSTRACT 

We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral 
resolution ^^CO (1-0) line survey of the inner ~ 10 x 6kpc of the M51 system, and the first wide-field imaging 
of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular 
Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer 
(PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1^' -resolution over the ~ 270^' x 
170^' field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, 
and brightness sensitivity of 0.4 K (la) in each 5 km s~^ -wide channel map. We find a CO-luminosity of 9 x 
10^Kkms~^pc^, corresponding to a molecular gas mass of 4 x 10^ M© for a standard CO-to-H2 conversion 
factor. Unexpectedly, we find that a large fraction, (50 ± 10)%, of this emission arises mostly from spatial 
scales larger than 3&^ c^ 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does 
not result from a processing artifact. We discuss its origin in light of the stellar component, the ^^CO/^^CO 
ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + 
IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a 
typical scale height of ~ 200 pc, substructured in unresolved filaments which fills ~ 0.1% of the volume. 

Keywords: galaxies: individual (M51) 



1. INTRODUCTION 

Along the path leading from the accretion of hot ionized gas 
onto galaxies to the birth of stars, the formation and evolution 
of Giant Molecular Clouds (GMCs) is the least well under- 
stood step. For example, the dependence of their mass distri- 
butions, lifetimes, and star formation efficiencies on galactic 
environment (e.g. arm, in terarm, nucle ar region) is largely un- 
known (for a review, see |McKee & O striker 2007). Because 
the Sun's position within the Milky Way disk makes GMC 
studies difficult within our own Galaxy, observations of GMC 
populations in nearby face-on galaxies offer the best way to 
address many of these unknowns. 

Complete CO maps that resolve individual GMCs have 
been carried out across the Local Group, allowing for the con- 
struction of mass functions and an estimation of G MC life- 
times via c omparison with maps at other wavelengths fKawa-j 
mura et al. 2009). To date, these observations have probea 
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mostly low m ass galaxies where Hi dominates the interstellar 
medium (g.g., [Blitz et al.|2007 ). The main reason is that the 
angular resolution re quired to identify i ndividual GMCs (typ- 
ical size ~ 40pc, g.g.,IS olomon et al.fT98 7) in any galaxy out- 
side the Local Group is extremely difficult to achieve. Reach- 
ing such resolutions with single dish telescopes remains im- 
possible in all but the very closest galaxies. This presents a 
major obstacle in linking our understanding of star formation 
and galactic evolution. Even for M31, the closest massive spi- 
ral galaxy to the Milky Way, the ^^CO (1-0) IRAM-30m map 
achieves a spatial resolution of only 8 5 pc and it suffers from 
projection effects ( [Nieten et al.||2006| . This is an important 
problem because these massive star forming spirals dominate 
the mass and light budget of blue galaxies and they h ost most 
of the star format ion in the present-day universe (e.g. Schimi- 
[novich et aLpQQT ). 

To remedy this situation, we used the Plateau de Bure Inter- 
ferometer (PdBI) to carry out the PdBI Arcsecond Whirlpool 
Survey (PAWS, Schinnerer e t al.|2013| ). The high quality re- 
ceivers and good weather conditions allowed PAWS to map 
the central, molecule-bright part of M51 at a resolution of 
~ 1.1^' (~ 40 pc) while still maintaining good brightness sen- 
sitivity (RMS ~ 0.4 K). M51 is the best target for such a pro- 
gram (see Table [T]). It is one of the closest (D ~ 7.6Mpc), 
face-on (/ '^ 21°) grand design spirals, and it has been exten- 
sively studied at essentially all wavelengths. In contrast to 
Local Group galaxies with a resolved GMC population, the 
molecular gas clearly domi nates the interstellar medium in- 
side the mapped region (g.g.,|Garcia-Burillo et al . 1993 ; Aalto 
'et al.'1999VSchuster et al. "2007 [ [Hitschf eld et al.|2009nLeroy 
et al. 2009; Koda et al. 2011). M51 thus offers the opportu- 
nity to relate the physical properties of molecular gas to spiral 
structure. We complemented the interferometric data with a 
sensitive (RMS ~ 16mK) map of the whole M51 system with 
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Figure 1. ^^CO (1-0) integrated emission of the inner ~ lOkpc x 6kpc of the galaxy NGC 5194 (M51a). The coordinate offsets are relative to the nucleus of 
NGC 5194 (see Table[T}. This image results from the joint deconvolution of the IRAM-30m single-dish and Plateau de Bure interferometer data sets. The image 
inserted at the bottom left is the ^^CO (1-0) integrated emission of the full M51 system (i.e., NGC 5194 + NGC 5195) as observed by the IRAM-30m telescope. 
The white horizontal rectangle shows the PAWS field of view. The images were scaled such that the angular resolution of both data sets occupies the same size 
on paper. In other words, the PAWS image shows the center of the small image zoomed by a factor of 21. 



Table 1 

Parameters for NGC 5194. 



Parameter 



NGC 5194 



Notes 



Ref. 



06| 



Morphological type 
Activity type 
Kinematic center 

Distance 

Systemic velocity 
Mean CO inclination 
Mean position angle 
Emitting surface 
Total CO luminosity^ 
Total molecular mass^ 
Mean brightness^ 
Mean mass surface density^ 



SA(s)bc pec 

Seyfert 2 

13h29m52 7087^+47° 11M2.789'' 

1.80'^0.81'' 

7.6±lMpc 

471.7±0.3kms-i 

21±3° 

173 ±3° 

1.9 X lO^pc^ 

1.4 X lO^Kkms-ipc^ 

6.2 X lO^M© 

7.6Kkms-i 

33M0pc^ 



Offset from phase center 

l'' = 37±5pc 
LSR, radio convention 



Wco > 3a 

in [LSR-120,LSR+120]kms- 

Helium included 



de Vaucouleurs et al. ( 199 
Veron-Cetty & Veron ( 2006| 
|HagiwaraH2007| 

CiarduUo et al. (2002} 
Shetty et al. (2007 1 
Colombo etal. (2013^ 
Colombo et al. ( 2013^ 
This work, Appendix A 
This work. Appendix A 
This work. Appendix A 
This work. Appendix A 
This work. Appendix A 



^ Using the IRAM-30m data. ^ The mean brightness and mass surface density are computed using the area with significant emission, i.e., Wqo > 3cr. 
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Figure 2. Pointing pattern of the two observed PdBI mosaics overlaid on the 
integrated emission of the ^^CO (1-0) line observed with the IRAM-30m. 
The 30 fields of each mosaic are displayed as red and blue circles, respec- 



or ~ 1430km s ^ at a spectral resolution of 1.25 MHz or 
3.25kms-i. 

We observed two 30-field mosaics, as described in Table [2] 
and shown in Fig. |2] Both mosaics were centered such that 
their combination covers the inner part of M51. The total 
field of view is approximately 270'' x 170''. Each pointing 
was observed during 3 x 15 seconds in turn. This allowed 
us 1) to observe one mosaic between two calibrations, which 
were taken every 22.5 minutes and 2) to minimize the dead- 
times due to moves from one field position to the next, while 
ensuring that the integration time per visibility (15 seconds) is 
short enough to avoid mixing independent uv plane informa- 
tion in all the configurations (s ee, e.g., Appendix C.l of ,Pety] 
|& Rodriguez-Femandez|201Ql for detailed calculations). An 
inconvenient aspect of such an observing strategy is that we 
obtained two data sets, observed in slightly different condi- 
tions, implying slightly different noise properties and uv cov- 
erage (i.e. slightly different dirty/synthesized beams). 

The field positions followed an hexagonal pattern, each 
field pointing being separated from its nearest neighbors by 
the primary beam full width at half maximum (FWHM), 
^fwhm- Nyquist sampling requires a distance between two 
consecutive pointings of X/dpnm along two orthogonal axes. 



wavelength and d^nm is the diameter of an interferometer antenna, Le 
for the PdBI. In our case, we have A/Jpnm ~ 36^^ 



tively. The diameter of each circle is A/Jprim, where A is the observation where A IS the observation wavelength and d^rim IS the diame 

15m terof the interferometer antennas. At PdBI, we typically have 

^fwhm = l-2A/(iprim. The hexagonal pattern used here thus en- 
sures Nyquist sampling along the Declination axis but a slight 
undersampling along the Right Ascension axis. 

The field-of-view was observed for about 169 hours of tele- 
scope time with 5 antennas in configuration D (19 hours) and 
6 antennas in configuration C (18 hours), B (57 hours) and A 
(75 hours). Taking into account the time for calibration and 
the data filtering applied, this translated into final on-source 
integration times (computed for a 6-antenna array) of useful 
data of 8.3 hours in D configuration, 15.2 hours in C config- 
uration, 43 hours in B configuration and 60 hours in A con- 
figuration. In each configuration, the time was approximately 
equally distributed between both mosaics. 



the IRAM-30m single-dish telescope. This allowed us to pro- 
duce a hybrid synthesis map — a joint deconvolution of the 
PdBI and IRAM-30m data sets — that is sensitive to all spa- 
tial scales between our synthesized beam and the PAWS field 
of view (see Fig.[T]). 

In Section |2] we detail the observing strategy and the data 
reduction. In Section [3j we show that a large portion (~ 
50% ± 10%) of the emission in our hybrid maps arises from 
faint, extended structures. We provide a detailed discussion of 
the nature of the gas responsible for this emission. We sum- 
marize our conclusions in Sect. |4] The Appendices provide 
details on technical aspects of our observations, reductions 
and analysis, and additional supplementary Tables and Fig- 
ures. 

2. PAWS DATA ACQUISITION AND REDUCTION 

This section presents the observin g stra teg y, the data reduc- 

focus on the 



tion and the resulting data set. Sects. |2.l| and [Z2 

tively 



23] explains 



PdBI and IRAM-30m data, respectively. Sect. 

how we combined these data to produce a final set of hybrid 

maps sensitive to all spatial scales. 

2.1. IRAM Plateau de Bure Interferometer data 

After a discussion of the observing setup, we describe the 
calibration of the interferometric data. 

2.1.1. Observations 

PdBI observations dedicated to this project were carried out 
with either 5 or 6 antennas in the A, B, C, and D configura- 
tions (baseline lengths from 24 m to 760 m) from August 2009 
to March 2010. The two polarizations of the single-sideband 
receivers were tuned at 115.090GHz, i.e., the ^^CO (1-0) 
rest frequency redshifted to the LSR velocity (471.7kms~^) 
of M51. Four correlator bands of 160 MHz per polariza- 
tion were concatenated to cover a bandwidth of ~ 550 MHz 



2.1.2. Calibration 

Standard calibration methods implemented inside the 
GILDAS/CLIC software were used for the PdBI data. The 
radio-frequency bandpass was calibrated using observations 
of two bright (-10 Jy) quasars, 0851+202 and 3C279, lead- 
ing to an excellent bandpass accuracy (phase rms ^1°, am- 
plitude rms ^1%). The temporal phase and amplitude gains 
were obtained from spline fits through regular measurements 
of the following close-by quasars: 1418+546, 1308+326, 
Jl 332+473. The flux scale was determined against the pri- 
mary flux calibrator, MWC349. The resulting fluxes of the 
calibration quasars are summarized in Table [TT] The absolute 
flux accuracy is r^ 10%. 

The data were filtered using statistical quality criteria on the 
pointing, flux, amplitude and phase calibrators. The source 
data were flagged when the surrounding calibrator measure- 
ments implied a phase rms larger than 40° , an amplitude loss 
larger than 22%, a pointing error larger than 30% of the pri- 
mary beamwidth and/or a focus error larger than 30% of the 
wavelength. Finally, the data were also flagged when the 
tracking error was larger than 10% of the field of view. This 
reduces the amount of usable data to 39%P^ 70%, 71%, and 

^^ This number takes into account visibilities that were flagged because of 
shadowing. 
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Table 2 

Parameters of the PdBI observations 














Molecule 


Transition Frequencies [GHz] 
Rest Tuned 


LSR 


Velocity [kms ^] 
Tuned Resolution 






i^CO 


(1-0) 115.271202 115.090 


471.7 







5 








Mosaic # 


Projection center (J2000) 
R.A Dec. 


A^fields 


Beam 




PA 

deg 








Top 

Bottom 

Full 


13^29'^54.92P 47° 12^21.589'' 
13^29'"50.504^ 47°11'03.642'' 
13^29'"52.532^ 47°llMl.982'' 


30 
30 
60 


1.12x0.92 
1.20 X 1.01 
1.16x0.97 


61.5 
84.0 
73.0 








Table 3 

Parameters of the 30m observations 










Molecule 
& Transition 


Frequency 
GHz 


Feff 5eff Resol.'' Resol.^ 
kms~^ arcsec 


Map Size 
arcmin^ 


Time'^ 
hours 


K[r;] 


(J 
mK [rmb] 


i^CO (1-0) 
i^CO (1-0) 


115.271202 
110.201354 


0.95 0.75 5.20/5.00 21.3/22.5 
0.95 0.76 5.44/5.00 22.3/23.6 


53 (- 
53 (~ 


6x10) 
6x10) 


17.3/41 
17.3/41 


285 
140 


16 

7.5 



^ The two values correspond to the backend natural channel spacing and to the channel spacing used to match the PdBI channel spacing. 

' The two values correspond to the natural FWHM of the beam and to the map resolution after gridding through convolution with a Gaussian. 

^ Two values are given for the integration time: the on-source time and the telescope time. 



Table 4 

Median noise levels of the different ^^CO (1-0) cubes. 



Cube type 




Resolution 










Icr Noise levels 














arcsec 


(mK[r„,b]) 




(Kkms-i) 






(M0PC-2) 












5kms-i 


25kms-i 


50kms-^ 


^ lOOkms-i 


5kms-i 


25kms-i 


50kms-i 


lOOkms-i 


IRAM-30m 




22.5 


16 


0.08 


0.18 


0.25 


0.36 


0.35 


0.78 


1.11 


1.57 


Hybrid synthesis 




6.0 


35 


0.17 


0.39 


0.55 


0.78 


0.76 


1.70 


2.40 


3.40 


Hybrid synthesis 




3.0 


106 


0.53 


1.18 


1.67 


2.37 


2.30 


5.15 


7.28 


10.30 


Hybrid synthesis 




1.1 


394 


1.97 


4.40 


6.23 


8.81 


8.56 


19.15 


27.08 


38.30 


PdBI-only 




6.0 


38 


0.19 


0.43 


0.60 


0.85 


0.83 


1.85 


2.62 


3.71 


PdBI-only 




3.0 


95 


0.47 


1.06 


1.50 


2.11 


2.06 


4.60 


6.50 


9.19 


PdBI-only 




1.1 


396 


1.98 


4.43 


6.27 


8.86 


8.62 


19.26 


27.24 


38.53 


Hybrid synthesis - PdBI 


-only 


6.0 


14 


0.07 


0.15 


0.22 


0.31 


0.30 


0.66 


0.94 


1.33 


Hybrid synthesis - PdBI 


-only 


3.0 


25 


0.12 


0.27 


0.39 


0.55 


0.53 


1.19 


1.69 


2.39 


Hybrid synthesis - PdBI 


-only 


1.1 


33 


0.17 


0.37 


0.52 


0.74 


0.72 


1.61 


2.27 


3.22 



71% of the data obtained in the D, C, B, and A configurations, 
respectively. 

2.2. IRAM-30m single-dish data 

A multipHcative interferometer filters out the low spatial 
frequencies, i.e., spatially extended emission. We thus ob- 
served M51 with the IRAM-30m single dish telescope on 
May 18-22, 2010 in order to recover the low spatial frequency 
("short- and zero-spacing") information filtered out by the 
PdBI. We describe here the observing strategy and the cali- 
bration, baselining and gridding methods we used to obtain 
single-dish data whose quality matches the interferometric 
data. 

2.2.1. Observations 

Table |3] summarizes the IRAM-30m observations. We used 
the EMIR receivers to map the ^^CO (1-0) and ^^CO (1-0) 
lines over a ~ 60 square arcminute field-of-view covering the 
M51 system, i.e., NGC5194 and its companion NGC5195. 
The upper sideband of the 3 mm separated sideband EMIR 
mixers (E090) was tuned at the ^^CO (1-0) frequency. The 
full 8 GHz bandwidth of the upper sideband was then con- 
nected to the WILMA autocorrelator backend. This allowed 
us to simultaneously measure the ^^CO and ^^CO lines (at 



115.271 and 110.201GHz, respectively). The backend chan- 
nel spacing is 2 MHz, which translates into a velocity chan- 
nel spacing of 5.4 and 5.2kms~^ at 110 and 115 GHz, respec- 
tively. 

We observed the galaxy in seven different patches. Four of 
these covered the central 400^' x 400^' part of the galaxy in 
different ways. Three additional patches extended the cover- 
age to include the ends of the spiral arms and the companion. 
Conditions during the observations varied from "good" sum- 
mer weather (~ 4 mm of precipitable water vapor) during the 
first three nights to "average" summer weather (r^ 7 mm of 
water vapor) over the last two nights. 

We used the position-switch on-the-fly observing mode, 
covering each field with back-and- forth scans along either the 
right ascension or declination axes. We slewed at a speed 
of ^ 8'7sec and we dumped data to disk every 0.5 seconds, 
yielding about 5.5 integrations per beam in the scanning di- 
rection (the HPBW of the IRAM-30m telescope at the fre- 
quency of the CO (1-0) line is ^ 22"). The scan legs were 
separated by 8'^ yielding Nyquist sampling transverse to the 
scan direction at 2.6mm. Each position in the central part was 
observed 34 times on average, with observations split evenly 
between right ascension- and declination-oriented scanning to 
suppress scan artifacts. The sky positions at the far end of the 
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Figure 3. Spatial distributions of the rms noise, peak intensity, integrated intensity and centroid velocity of the ^^CO (upper panels) and ^^CO (lower panels) 
(1-0) emission as observed with the IRAM-30m telescope. The angular resolution is displayed as a circle in the bottom left corner of each panel. The intensity 
scale is shown on the right side of each panel with units indicated in the caption title. The major and minor kinematic axes are displayed as perpendicular dotted 
lines. The dotted circles show the two inner corotation resonances at radii equal to Th" and 55^^ while the dashed circle marks the start of the material arms at a 
radius equal to 85^^ ^Meidt et al.|2013^ . 



CO spiral arms were observed 12 times so that the final ef- 
fective integration time on the extensions is somewhat shorter 
than on the main field (see Fig. [3]). 

We observed the hot and cold loads plus the sky contribu- 
tion every 12 minutes to establish the temperature scale and 
checked the pointing and focus every ~ 1 and r^ 4 hours. The 
IRAM-30m position accuracy is ~ 1" . 

2.2.2. Calibration and gridding 

We reduced the IRAM-30 m data usin g a combination of the 
G I LDA^ software suite (Pety"2005') and an IDL pi peline 
developea for the IRAM HERACLES Large Program (Leroy 
let al.|2QQ9D . 

First, we calibrated the temperature scale of the data in 
GILDAS/ MIRA based on the hot an d cold loads plus sky ob- 
servations ( [Penzias & B urrus||1 973|). The resulti ng flux ac- 
curacy is better than 10% (Kra mer et aTlpOOSI ). We then 
subtracted the "OFF" spectrum from each on-source spec- 
trum. We used GILDAS/ CLASS to write these calibrated, 
off- subtracted spectra to F I T S tables which we read into IDL 
for further processing. 

In IDL, visual inspection indicated the presence of signal 
in the [-200, 3 00 km s~^] velocity range around the systemic 
velocity of the galaxy. About 1/1 6th of the 8 GHz bandwidth 
(i.e. about 1000 km s~^) centered on the line rest frequency 
were thus extracted from the calibrated spectra to reduce the 
computation load of the next data reduction steps. We then 
fit and subtracted a third-order baseline from each spectrum. 



When conducting these fits, we use an outlier-resistant ap- 
proach and exclude regions of the spectrum that we know to 
contain bright emission based on previous reduction or other 
observations. We experimented with higher and lower order 
baselines and found a third-degree fit to yield the best results. 
After fitting, we compared the rms noise about the baseline fit 
in signal-free regions of each spectrum to the expected theo- 
retical noise and used this to reject a few pathological spectra. 
For the most part the data are very well behaved and this is a 
minor step. 

We gridded the calibrated, off-subtracted, baseline- 
subtracted spectra into a data cube whose pixel size is 4'^ 
Doing so, we weighted each spectrum by the inverse of the 
associated rms noise. We employed a gaussian convolution 
(gridding) kern el with FWHM 7.V \ ~ 1/3 of the primary 
beam (e.g., see |Mangum et al.|2007| ). This gridding increases 
the FWHM of the effective beam from - 21.3'' to - 22.5'' at 
115 GHz. 

After gridding, we fit a second set of third-order polyno- 
mial baselines to each line of sight through the cube. The pro- 
cess of the initial fitting and gridding is linear, so that these 
fits represent refinements to our initial fits after the averag- 
ing involved in gridding. We experimented with seve ral more 
advanced pr ocessing options such as PLAITing (Em erson & 
Graeve"1988) and flagging of standing waves. However the 



^^ See http : / /www, ira m. f r/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 



data were very clean and none of these algorithms improved 
the quality. 

Figure [3] presents the reduced, calibrated, gridded IRAM- 
30m maps of the ^^CO (1-0) and ^^CO (1-0) line emis- 
sion. The figure displays the spatial distributions of noise. 
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Figure 4. Top, left: Area of the CLEAN mask, normalized to its maximum 
value, as a function of the channel velocity. Top, right: Ellipses representing 
the mean Gaussian beams for the top (red) and bottom (blue) mosaics. The 
black ellipse represents the Gaussian beam used for the restoration of the 
CLEAN component list. Middle: uv coverages for the fields of the top (left 
panel) and bottom (right panel) mosaics. Bottom: Zoom of the dirty beams 
for the top (left panel) and bottom (right panel) mosaics. 



peak temperature, integrated emission and centroid velocity. 
Fig. [26l and [27] (available in the electronic version only) dis- 
play tEe channel maps of these lines. 



2.3. Combination, imaging and de convolution 

The interferometric and single-dish data provide us with 
two data sets, which sample the high and low spatial frequen- 
cies, respectively. It is thus possible to produce two different 
deconvolved results: 1) one obtained from the interferomet- 
ric data set alone, and 2) one obtained from the combination 
of the interferometric and single-dish data. While the latter 
is the desired final product, as it is sensitive to all measured 
spatial scales, the former is often produced because no single- 
dish measurements will be available or because they have not 
yet been acquired. In this paper, we will present both data 
cubes to emphasize the amount of flux which is recovered in 
the interferometric-only data set. Moreover, the angular res- 
olution of the interferometric data is not uniquely defined. It 
depends on the weighting scheme chosen. For instance, it 
is sometimes useful to produce data cubes at lower angular 
resolutions to improve the brightness sensitivity, i.e., the sen- 
sitivity to extended emission. We exploit this here in addition 
to producing the full resolution cube. 

This section explains 1) the generic imaging and deconvo- 
lution methods used to produce all these data cubes, 2) how 
these methods influence the amount of flux recovered in the 
interferometric-only data, and 3) the additional steps required 



to image jointly the single-dish and interferometric data. 

2.3.1. Generic methods 

Each interferometric pointing was imaged and a single dirty 
image was built by linear combination of the 60 individual 
dirty images. The dirty image was then deconvolved us- 
ing an adaption of the Hogbom CLEAN algorithm. A de- 
tailed account of the GILDAS/MAPPING implementation of 
the imag ing and deconvolution processing of m osaics can be 
found in |Pety & Rodriguez-Fernandez] pOlO] ). To help the 
deconvolution, masks indicating the region where to search 
for CLEAN components were defined on individual channels 
from the short-spacing data cube. This cube was convolved 
with a Gaussian kernel to a final angular resolution of 30'^ 
The CLEAN masks were then defined by all the pixels whose 
signal-to-noise ratio was larger than 6. This method was de- 
signed to avoid biasing the deconvolution by defining masks 
wide enough to encompass all the detected signal from M51. 

The deconvolution of each channel was stopped when a 
fraction of the maximum number of clean components were 
found. This fraction was defined as the ratio of the area of the 
current channel mask to the area of the wider channel mask 
(see left panel of Fig.ffl. The deconvolution was assumed to 
have converged under three conditions. First, the cumulative 
flux as a function of the number of clean components con- 
verged in each channel. Second, the residual channel images 
looked like noise. Both criteria indicated a satisfying conver- 
gence of the deconvolution. Finally, we deconvolved again 
the data using exactly the same method except that we dou- 
bled the maximum number of clean components. The sub- 
traction of both cubes looks like noise. 

The list of clean components was regularized with a Gaus- 
sian beam and the residual image was added to obtain the final 
cube. As all the 30 fields of the top and bottom mosaics were 
regularly observed in short cycles of 22.5 minutes, the synthe- 
sized beams do not vary inside each mosaic. However, the two 
mosaics were observed at different times implying a slight dif- 
ference of the synthesized beam between both mosaics (see 
Table [2]). We used the same averaged Gaussian restoration 
beam ror the northern and southern mosaic. This process is 
valid because 1) the Gaussian fits of the synthesized beams 
are similar as shown in the top right panel of Fig. |4| and 2) 
the remaining flux in the residual image is negligible or unde- 
tectable {i.e., below the noise limit). The resulting data cube 
was then scaled from Jy/beam to Tiib temperature scale using 
the restoration beam size. 

While the natural velocity channel spacing of the inter- 
ferometer backend is 3.25 km s~^ at 2.6mm, we smoothed 
the data to a velocity resolution of 5kms~^ This decreases 
the effect of correlation between adjacent frequency channels 
output by the correlator. This also increases the signal-to- 
noise ratio per channel (an important factor for the deconvo- 
lution) and speeds up the processing. Signal is present be- 
tween -110 and +110kms"\ relative to NGC5194's LSR 
systemic velocity of 47 1.7 km s~^ We thus imaged and 
deconvolved 120 channels producing a velocity range of 
[-297.5, +297.5 km s~^], implying that about two-thirds of the 
channels are devoid of signal. On a machine with 2 octa- 
core processors and 72 GB of total RAM memory, the de- 
convolution of the 120 channels up to a maximum number of 
320000 clean components typically took 38.5 hours (human 
time). The deconvolution duration increases linearly with the 
maximum number of clean components per channel. 
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Figure 5. Deconvolved flux as a function of velocity for the hybrid synthesis 
(top) and the PdBI only (bottom) data sets. Yellow, red, and black curves 
correspond to cubes imaged at an angular resolution of 1" , 3", and 6", re- 
spectively. The percentage of total flux recovered compared to the 6" cube is 
indicated in the top right corner using the same color coding. 



2.3.2. PdBI only 

To start, we imaged and deconvolved the PdBI data with- 
out the short- spacings from the IRAM-30m observations. 
Achieving the convergence of the deconvolution algorithm at 
a given angular resolution is insufficient to prove that all the 
flux was recovered in PdBI-only data sets. Indeed, the ab- 
sence of zero spacing implies that the total flux of the dirty 
image is zero valued, and it is the deconvolution algorithm 
which tries to recover the correct flux of the source. This 
works only when the source is small compared to the primary 
beam and the signal-to-noise is large enough. No deconvo- 
lution algorithm will succeed in recovering the exact flux of 
even a point source at low signal-to-noise. Indeed, the algo- 
rithm recovers only the flux which is above a few times the 
noise rms. Adding the deconvolution residuals will not help 
because the dirty beam (i.e., the interferometer response) has 
a zero valued integral, i.e., the residuals always contain zero 
flux. 

Moreover, a given interferometer needs 2"^ as much observ- 
ing time to keep the same brightness sensitivity when just 
doubling the angular resolution (assuming similar observing 
conditions). Such an increase in observing time is impracti- 
cal. The brightness sensitivity thus decreases quickly when 
the angular resolution improves. In the PAWS case, we ap- 
proximately doubled the observing time every time we went 
to the next wider interferometer configuration, which typi- 
cally doubled the angular resolution. This allowed us to reach 
a median noise of 0.4 K at full resolution, i.e. 1 . 16^' x 0.97^' at 
a position angle (PA) of 73 deg (when using natural weighting 
of the visibilities). While this is the best (pre- ALMA) sensivi- 
tivity reachable for such a large mosaic, it is also much higher 
than the sensitivity of 16mK we reached with the IRAM-30m 
at an angular resolution of 22.5'^ This may mean that faint in- 
tensities may be hidden in the noise of the 1^' -resolution cube. 

Multi-resolution CLEAN algorithms (only available for 



single-field observations in GILD AS and thus not used for the 
PAWS mosaic) partly solve this problem of brightness signal- 
to-noise ratio because the deconvolution simultaneously hap- 
pens on dirty images at different resolutions and thus differ- 
ent brightness noise levels. Interferometric brightness noise 
is a compromise between the synthesized angular resolution 
and the time spent in the different configurations. To check 
what happens with our deconvolution algorithm, we tapered 
the visibility weights to increase the brightness sensitivity at 
the cost of losing angular resolution, as this method is to first 
order similar to a Gaussian smoothing in the image plane. 
We choose Gaussian tapering functions such that we obtained 
synthesized resolutions of 3^' and 6'^ TableHlists the typical 
rms noise values for the 3 different resolutions, i.e., 0.4, 0.1 
and 0.03 K at respectively V\ ?>", and 6". 

The lower panel of Fig. [5| visualizes the flux found by the 
CLEAN algorithm at the 3 aifferent resolutions as a function 
of velocity. The Hogbom CLEAN algorithm finds 40% more 
flux at y than dXl" . On the other hand, only 4% more flux is 
recovered at 6" than at ?>" even though the brightness sensitiv- 
ity increases by a factor of ~ 2.5. Since the typical resolution 
of PdBI at 3 mm in its most compact configuration is 6" , it 
means that PdBI reaches its maximum brightness sensitivity 
in this configuration. Hence, further tapering the data would 
not recover more flux. Recovering only a marginal additional 
amount of flux when going from ?>" to 6" thus implies that we 
recovered at these resolution all the flux present in the inter- 
ferometric data. 

2.3.3. Hybrid synthesis (PdBI + IRAM-SOm) 

The hybrid synthesis is a joint deconvolution of the PdBI 
and IRAM-30m data sets. The IRAM-30m and PdBI data 
sets were first made consistent using the following 4 steps. 1) 
We converted the IRAM-30m spectra to main beam tempera- 
tures (r^b) using the forward and main beam efficiencies (Feff 
& ^eff) given in Table [3] 2) These spectra were reprojected 
on the projection center used for the interferometric data set. 
3) The LSR systemic velocity was set to zero to mimic ob- 
servations at the redshifted frequency. 4) The velocity axis 
was resampled to 5kms~^ channel spacing. The last transfor- 
mation introduced some correlation between channels as the 
autocorrelator natural channel spacing is only 5.2kms~^ at the 
^^CO (1-0) fr equency. 
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GILDAS/MAPPING software and 
from the IRAM-30m were used to create the short-spacing 
visibilities not sampled by the Plateau de Bure interferometer. 
In short, the maps were deconvolved from the IRAM-30m 
beam in the Fourier plane before multiplication by the PdBI 
primary beam in the image plane. After a last Fourier 
transform, pseudo-visibilities were sampled between and 
15 m (the diameter of the PdBI antenna). These visibilities 
were then merged with the interferometric observations. The 
relative weight of the single-dish versus interferometric data 
was computed in order to get a combined weight density in 
the uv plane close to Gaussian. Since the Fourier transform 
of a Gaussian is a Gaussian, and the dirty beam is the Fourier 
transform of the weight density, this ensures that the dirty 
beam is as close as possible to a Gaussian, making this 
criterion optimum from the deconvolution point-of-view. In 
general, the short spacing frequencies are small compared to 
the largest spatial frequency measured by an interferometer. 
This implies we can use the linear approximation of a 
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Figure 6. Spectra, averaged over the field of view, of the hybrid synthesis 
(blue), PdBI-only (red) and filtered (= hybrid-PdBI, white) emission. The 
percentage of flux recovered in each spectrum is written in the top right cor- 
ners with the same color code. The computations were done at the angular 
resolution displayed at the top left corner of each panel. 

Gaussian in the vicinity of its maximum, i.e., one can assume 
the Gaussian is constant in the range of frequencies used for 
the processing of the short spacings. We thus need to match 
the single-dish and interferometric densities of weights in the 
inner region of the uv plane. 

In practice, we compute the density of weights from the 
single-dish in a uv circle of radius 1.25Jprini and we match it 
to the averaged density of weights from the interferometer in 
a uv ring between 1.25 and 2.5Jprim. Experience shows that 
this gives the right order of magnitude for the relative weight 
and that a large range of relative weight aroun d this value 
gives very similar final results (see e.g. Fig. 5 of [Rodriguez- 1 
[Fernande z et al.|[20Q8| ). In our case, this computation was 
independently done for each interferometric pointing to take 
into account the fact that the two mosaics used to produce 
the final image were observed in slightly different conditions. 
The relative weight varies typically by 1% from pointing to 
pointing in each mosaic, and by ~ 6-7% between both mo- 
saics. 

Contrary to interferometric-only data sets, the upper panel 
of Fig. [5] shows that the deconvolved flux in a hybrid synthe- 
sis is the same for the 3 different angular resolutions. In other 
words, the deconvolved flux is independent of the brightness 
sensitivity reached, as it is fully constrained by the zero spac- 
ing amplitude. Having high signal-to-noise zero spacing data 
thus ensures that the total flux inside the deconvolved cube 



will be the total fl ux of t he single-dish data in the same field 
of view. Section 13.1.31 checks this for the PAWS data set. 
This is linked to the fact that the dirty beam integral is now 
normalized to unity. Both the dirty image and the residual 
fluxes are meaningful in this case. This enables an additional 
check of the convergence of the deconvolution algorithm for 
the hybrid synthesis deconvolution: We checked that less than 
1.2% of the clean flux remains in the residual cube in the 
[-100,+100]kms~^ velocity range, even at \" angular reso- 
lution. 

Fig. [28] and [29] (available in the electronic version 
only) display the channel maps of the 1^' -resolution hybrid 
(30m+PdBI) synthesis cube at negative and positive veloci- 
ties (relative to the systemic velocity), respectively. 

3. A LUMINOUS COMPONENT OF EXTENDED 
EMISSION 

This section presents the unexpected finding that about half 
the CO luminosity in the hybrid map ar ises from a faint, ex- 
tended component. Sect. [3.1 1 and |3.1.5| demonstrate that this 
result is unlikely to reflect an artifact m the data. Sect.|3.2|ex- 



plores how the emission in the hybrid map breaks apart into a 
bright, compact and a faint, extended component. It al so com- 
pares the structures of these distributions. Sect. 3.3 and |3.4| 
consider the nature of the CO emission components m light of 
1) calculations of verti cal d isk structure, and 2) the observed 
^^CO/^^CO ratio. Sect. 3^ discusses our interpretation of this 
emission component. 

3.1. Verifying the existence of an extended component 

This section shows that (50 ± 10)% of the total flux in the 
PAWS field of view is filtered out by the interferometer. After 
a thorough discussion of several checks to prove that this is 
a true effect, we show that the filtered emission has typical 
spatial scales larger than 36'^ i.e., 1.3kpc(Z)/7.6Mpc), where 
D is the assumed distance of M5 1 . 

3.1.1. Flux recovered in the PdBI-only data cubes 

Fig. J6] compares the spectra averaged over the field of view 
of 1) the hybrid synthesis data set (in blue) and 2) the PdBI- 
only data set (in red) for the 1, 3, and 6" data cubes from top 
to bottom. These spectra were obtained as a mean over the 
PAWS field of view contracted by about 1 primary beam width 
to decrease the influence of increasing noise at the map edges. 
The differences between the hybrid synthesis and PdBI-only 
spectra are displayed in white. 

The main result is that only half of the total flux is recov- 
ered in the PdBI-only data set. Indeed, only 37% of the to- 
tal flux is recovered at \" but this is attributed to the rela- 
tively "low" brightness sensitivity reached at this resolution 
(see Section |2.3.2| ). Approximately 50% of the total flux is 
recovered both at ?>" and 6" , while the brightness sensitivity 
differs by a factor of 2.5 between them. We stress that al- 
though the deconvolution recovers 44% more flux at 6" than 
at \" , the difference amounts to only 13% of the total flux 
present in the PAWS field of view. 



3.1.2. Coordinate registration 

We checked the overall registration of the IRAM-30m 
moment-0 images against the PdBI-only data. We convolved 
these "reference" maps to the resolution of the IRAM-30m 
data. We then repeatedly shifted the images relative to one an- 
other and recorded the cross-correlation between the IRAM- 
30m data and the other images. The overall registration of 
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Figure 7. Channel by channel flux comparison of the PdBI+30m, IRAM- 
30m, CARMA and BIMA data cubes over the PAWS field of view. The 
bottom panel shows the mean CO temperature as a function of the velocity 
and the top panel shows the relative difference with respect to the PdBI+30m 
data set. The dotted horizontal line indicates perfect agreement. 



Table 5 

Comparison of registration and total flux within different ^^CO (1-0) 
datasets. 



Data set 


Offset 


Integrated flux 


Flux difference 




arcsec 


lO^Kkms- 


ipc2 


%ofPdBI+30mflux 


IRAM-30m 





7.82 




-0.5 


PdBI+30m 


< 1 


7.86 




+0.0 


CARMA+45m 


< 1 


7.41 




-5.7 


BIMA+12m 


< 1 


7.03 




-10.6 



the IRAM-30m with the PdBI data appears to agree within 
^ v. In additio n we obtained the same agreement with 
the B IMA SoNG ( [Heifer et al.| [2Q03) and CARMA (,Koda| 
et al.pOH ) data. This good agreement is expected because 



radio-observatories check the coordinate registration directly 
against the radio-quasars used to define the equatorial coordi- 
nate frame. 

3.1.3. Flux calibration 

As a test of our calibration strategy, we compared the flux 
within the hybrid synth esis data cube w ith recent CO surveys 
of M51 by BIMA (Hei fer etal.|[2003] ) and CARMA (Koda 
[etal.p Qll). Although conceptually straightforward, this com- 
parison must be done with care since each dataset has dif- 
ferent angular resolution, channel width, noise characteristics 
and field of view. To obtain the most meaningful comparison, 
we smoothed the BIMA, CARMA and hybrid synthesis data 
cubes to the same angular resolution as the IRAM-30m data 
(22f'5) and interpolated them onto a common (/,m,v) grid 
with square pixels of 4" size, and a velocity channel spacing 
of 5kms~^ The grid uses a global sinusoid (GLS) projection, 
and is centered at RA 13:29:54.09, Dec +47:11:38.0 (J2000). 

To define the region for the flux comparison, we con- 
structed a mask of significant emission within the IRAM-30m 



data cube using the dilated mask technique described in Ap- 
pendix |B] In summary, this mask is obtained by the contigu- 
ous extension down to an intensity level of 2 times the rms 
noise a of all the pixels whose intensity is above Scr. The 
IRAM-30m a's were estimated for each sightline using the 
median absolute deviation. The resulting mask was applied to 
all four data cubes. The total integrated flux within a spatial 
region corresponding to the PAWS field of view contracted by 
~ 1 primary beam width was then computed for each data set. 

The results are listed in Table [5l The flux of the combined 
hybrid synthesis cube agrees with the flux of the IRAM-30m 
to within 1%, which is expected since the algorithm that we 
have used to combine the interferometer and single-dish data 
is designed to conserve flux. The integrated fluxes estimated 
from the PdBI+30m, CARMA+45m and BIMA+12m surveys 
agree to within 10%, which is the typical accuracy of the abso- 
lute flux calibration of radio observatories at millimeter wave- 
lengths. 

Figure [7] shows the global spectrum of the four surveys 
from the comparison region (bottom panel) and the relative 
flux difference per velocity channel between PAWS and the 
three other data sets (top panel). In general, the flux differ- 
ences for individual velocity channels are consistent with the 
integrated measurements, i.e., there is ~10% (5%) less flux in 
each channel of the BIMA+12m (CARMA+45m) data cube 
compared to PdBI+30m. 

These tests confirm the IRAM-30m flux calibration accu- 
racy and its correct transmission to the hybrid synthesis data 
set. As a test of the flux accuracy of the PdBI-only data set, 
we compared the efficiencies resulting from the interferomet- 
ric flux calibration to the well-known PdBI antenna efficien- 
cies as measured through regular holographies. These agree 
within ~ 10%. The key point here is that the flux calibration 
is relatively straightforward at 3 mm, because the antennas 
have an excellent efficiency and the atmosphere was mostly 
transparent when the data was acquired (as confirmed by the 
achieved system temperatures listed in Table [10 ). 



3.1.4. Amplitude of the brightness Fourier transform as a function 
of the uv distance 

It is important to test how well the hybrid synthesis and the 
PdBI-only data cubes agree with each other at high spatial fre- 
quencies. To do this, Fig.[8]compares the azimuthal average of 
the amplitude of the Fourier transform of the PdBI-only and 
the hybrid synthesis data cubes at &^ resolution as a function 
of the uv distance for every third velocity channel between 
[-102.5, +107.5] kms~^ In such a representation, the ampli- 
tude at the zero uv radius is equal to the total flux present 
in the channel map. We thus normalized all the curves in- 
side each panel so that their zero- spacing values represent the 
fraction of the total flux available in the hybrid synthesis cube 
at this velocity channel. For reference, 1) the power spec- 
trum would be computed as the square of these curves, 2) 
the dashed vertical lines indicate the minimum radius directly 
measured by the interferometer, i.e., rmm — 15.1m. More- 
over, noise and signal behave differently after undergoing a 
Fourier transformation. Noise (e.g., signal-free channels at 
the edges of the velocity range in the figure) will give a shape 
consistent with the Fourier transform of the Gaussian restora- 
tion beam (see Chapter 17 of Bracewell 2000). On the other 
hand, signal channels display the product of Fourier transform 
of the true source brightness and the Fourier transform of the 
Gaussian restoration beam as a function of the uv spacing. 
Hence, channels with significant line signal fall more quickly 



10 



Pety et al. 



- PdBI-only data 

- Hybrid synthesis data 



Hybrid synthesis beam 

-87.5 -- 




UV radius [m] 



Figure 8. Variations of the azimuthal average of the Fourier transform am- 
plitude as a function of uv radius for every third velocity channel between 
-102.5 and +107.5 km s~^ . On the y-axis, we plot the percentage of the maxi- 
mum flux at the zero uv radius. The pink dashed curve displays the theoretical 
shape for a signal-free channel, i.e., the Fourier Transform of the hybrid syn- 
thesis beam. The pink and blue solid curves show, respectively, the hybrid 
synthesis and PdBI-only data imaged at 6'' resolution. The dotted vertical 
line indicates the minimum uv radius measured by the interferometer (i.e. 
~ 15 m). The dotted horizontal line represents a 50% reduction in flux. 
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Figure 9. Zoom of the variations of the azimuthal average of the Fourier 
transform amplitude as a function of uv radius for every third velocity channel 
between -102.5 and -i-107.5kms~^. The y-axis shows the percentage of the 
maximum flux at zero uv radius. The dashed pink and green curves display 
the Fourier Transforms of the hybrid synthesis and the IRAM-30m beams, 
respectively. The solid pink curve shows the hybrid synthesis data imaged at 
6" resolution. The solid green and cyan curves represent, respectively, the 
IRAM-30m data before and after deconvolution from the 22.5^^ 30m beam 
and convolution with a 6" beam. The dashed vertical line indicates the mini- 
mum uv radius measured by the interferometer (i.e., ~ 15m). 

than the Gaussian restoration beam. 

The good agreement between the ampHtudes from the hy- 
brid and PdBI-only data at radii larger than rmin provides con- 
fidence in the deconvolution results of the compact sources 
(whose angular extent is smaller than A/rmin ^ 35^0 in the 



PdBI-only data set even though the short- spacings are miss- 
ing. Below Tmin, the two curves diverge because the decon- 
volution can only "extrapolate" up to the zero uv radius for 
the PdBI-only data set, while it "interpolates" for the hybrid 
synthesis data set. As the amplitude at zero spacing is propor- 
tional to the total flux, these plots clearly indicate that PdBI- 
only data miss a large fraction of the flux even for the 6" res- 
olution cube, whose deconvolution is immune to low signal- 
to-noise effects. 

Moreover, the change of slope visible in Fig. [8] for the hy- 
brid synthesis amplitude between ~ 5 and 15 m is not a pro- 
cessing artifact. Indeed, Fig. [9] displays a zoom of the az- 
imuthal average of the amplitude of the Fourier Transform of 
the IRAM-30m and the hybrid synthesis data cubes at 22. 5'^ 
and &^ resolution (green and pink curves, respectively). These 
curves cannot be directly compared because they result from 
two different measurement equations (the single-dish and the 
interferometric ones). The solution is well-known: We just 
have to apply to the IRAM-30m curve the standard process- 
ing s teps n eeded to produce short- spacing visibilities (see 
Sect. |2.3.3| ). Indeed, after deconvolution of the single-dish 
brightness by the 22. 5^' single-dish beam and multiplication 
with the 6" interferometric primary beam, the amplitudes of 
the Fourier transform of both data sets agree well. Since a 
Gaussian of 22.5'' FWHM has a flux equal to 77, 45, and 21% 
of its maximum at 5, 10, and 15 m, the observed change of 
slope of the deconvolved 30m data cannot be attributed to im- 
precision in the deconvolution of the IRAM-30m beam. This 
change of slope in flux is responsible for the failure of the de- 
convolution to extrapolate the correct total flux at the zero uv 
radius for the PdBI-only data set. 

3.1.5. Impact of the IRAM-SOm error beam 

The beam efficiency of the IRAM-30m telescope at the fre- 
quency of the ^^CO (1-0) line is predicted to be 0.78, based 
on the interpolation of the measured efficiencies at 86, 145, 
210, 260, and 340 GHz using the Ruze formul^ Moreover, 
the forward efficiency is 0.95 at the same frequency. This 
implies that about 22% of the measured flux in a given di- 
rection results from beam pick-up from solid angles outside 
the main beam. Of this amount, 5% originates from the rear 
lobes, which mainly collects diffuse emission from the tele- 
scope backside. Hence, about 17% of the measured flux is 
picked up by the rings of the telescope diffraction pattern and 
the error beams due to the limited surface accuracy of the tele- 
scope. In Appendix JC) we assess in-depth how much these 
effects contribute to the extended emission. Here, we summa- 
rize the main results. 

We estimate that the peak brightness due to the error beam 
contribution is at most 55 mK, i.e., about 3.5 times the me- 
dian noise level of the IRAM-30m observations. The median 
value of the error beam contribution is 19mK, while pixels 
brighter than 5 times the noise level have a media n of 14 0mK 
in the extended emission measured at 6" (see Sect. 3.2.1 ). The 
emission associated with the error beams has a very differ- 
ent signature from that of extended emission both in space 
and velocity. In particular, the putative error beam contribu- 
tion translates into much wider lines than actually observed. 
Moreover, the typical angular scales of the error beams are 
large, implying that the flux is scattered at low brightness level 
over wide regions of the sky. Finally, we estimate that at most 



For details, see http://www.iram.es/IRAMES/mainWiki/ 
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Figure 10. Same as Fig.jslbut the pink dot-dashed curve displays the average 
noise level as a function or the uv radius. The red curve shows the difference 
between the hybrid synthesis and the PdBI-only at 6^' resolution. The 3 dot- 
ted vertical lines indicate dpnm, 2<iprim. and SJprim, where dp^m is the primary 
beam diameter of the interferometer. 

10% of the total flux in the PAWS field of view can be due to 
error beam scattering. 

Taken together, these arguments lead us to the conclusion 
that the IRAM-30m error beam cannot explain the presence 
and properties of the extended emission. 

3.1.6. About half the CO luminosity arises from an extended 
emission component 

We now consider why (50 ± 10)% of the total flux is missing 
in the PdBI-only data set. The direct interpretation is that this 
missing flux has been filtered out because the PdBI does not 
measure visibilities at spacings shorter than ~ 15m, i.e. the 
well-known short-spacing problem. We now wish to quantify 
the spatial scales at which the emission filtered by the PdBI is 
structured. To do this. Fig. [T0| compares as a function of the 
uv distance 1) the amplitude of the difference of the Fourier 
transforms between the hybrid synthesis and the PdBI-only 
data cubes (red plain curve), with 2) the noise of the hybrid 
synthesis Fourier transform (pink dashed curve), computed 
by averaging the signal-free channels. For all the velocity 
channels showing some signal, the amplitude of PdBI-filtered 
emission is (much) larger than the noise level only up to a 
radius of ~ 15 m. This implies that a very large fraction of 
the missing flux is structured only at spatial scales larger than 
~ A/rmin = 36'^ However, the difference and noise curves stay 
close to each others up to 35 m for a few velocity channels 
{e.g., V = -42.5 or +32.5 kms~^). This indicates there is prob- 
ably a small fraction of the missing flux which is structured 
at spatial scales between 15^' and 36'^ In the following, we 
will thus state that the missing flux is structured mostly at spa- 
tial scales larger than ~ A/15 m = 36'^ or 1 .3 kpc (Z)/7.6Mpc), 
where D is the distance of M5 1 . Conversely, this also means 
that the flux recovered in the PdBI-only cubes is structured 
mostly at spatial scales smaller than 36'^ 

3.2. Structure of compact and extended emissions 

In order to investigate the nature of the emission filtered by 
the PdBI interferometer, we need to image it. To do this, we 



subtracted at each angular resolution (1,3 and 6") the PdBI- 
only from the hybrid synthesis data cubes to determine the 
properties of the emission filtered by the PdBI. We were care- 
ful to image both data sets on the same spatial and spectral 
grid. We also used the same weighting scheme, deconvolu- 
tion method, stopping criterion and restoration beam. 

In the following, we will refer to the 3 cubes as hybrid 
synthesis (PdBI+30m), PdBI-only, and subtracted {i.e., hy- 
brid synthesis minus PdBI-only) cubes. The PdBI-only and 
subtracted cubes are equivalent to a decomposition of the hy- 
brid synthesis signal into two kinds of source morphology: 
Compact (angular scales ^36^0 and extended (angular scales 
>?>6") sources. This decomposition is a convenient instru- 
mental side effect. As such, it is arbitrary and only a multi- 
scale analysis would deliver the full spatial distribution of 
the emission. However, Fig. [8] and [9] show that the Fourier 
transform amplitude changes its slopebetween ~ 5 and 15 m. 
Thus, the intrinsic distribution of spatial scales of M51 is such 
that our decomposition makes sense. 

For completeness, we carried out this decomposition for the 
1, 3 and 6" cubes. While this gives a taste of multi-scale ap- 
proaches, two of these cubes stand out. The PdBI-only I" data 
cube is the one where the "compact" sources are best resolved 
while the subtracted 6" data cube is the one which gives the 
most accurate description of the extended component. Finally, 
we smoothed the PdBI-only 6" cube to the resolution of the 
IRAM-30m data cube (22.5^0 and we subtracted them from 
the 30m data. This allows us to obtain an idea of the relative 
contribution of both emission components at the typic al reso - 



lution of single-dish observations as discussed in Sect. 3.4.3 

In the rest of this section, we first describe the statistical 
distributions of noise and signal for the different data cubes, 
as this allows us to discuss how beam dilutions and signal-to- 
noise ratios evolve with angular scale. We then comment on 
the 2-dimensional spatial distributions of the line moments. 
We also present their azimuthal averages, which we will later 
use to analyse the vertical structure of both components. Fi- 
nally, we show the kinematics of the two components along 
the M51 major axis. 

3.2.1. Noise and signal distributions 

In order to quantify the (/,m,v)- volume filling factors and 
the beam dilution properties of the compact and extended 
emission, we now describe the distribution of the noise and 
brightness values. 

Figure 1 1 displays the cumulative histograms of the rms 
noise for me 4 different resolutions (1, 3, 6, and 22.5^0. and 
for the 3 kind of cubes: Hybrid synthesis (top), PdBI-only 
(middle) and subtracted (bottom). The left column displays 
the raw histograms while the right column presents the his- 
tograms normalized by their median value, as this eases the 
comparison of shapes. The histograms show that the noise 
distributions are well centered on the median value, implying 
relatively uniform noise properties. The IRAM-30m single- 
dish data display an increase of the histogram population at 
high noise values because the noise distributions are homo- 
geneous over the inner 400^' x 400^' field of view but they 
increase on the northern and southern patch. A similar effect 
is seen for the interferometric data, which is directly linked to 
a noise increase at the edges of the field because of the cor- 
rection of primary beam attenuation. Moreover, in the single- 
dish data, the ^^CO median noise level (16mK) is about twice 
as high as the ^^CO one (7.5 mK) because the ^^CO (1-0) line 
is closer to an oxygen atmospheric line. 
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Figure 11. Cumulative histograms of the rms noise for the ^^CO (1-0) 
(dashed line) and ^^CO (1-0) (solid lines) cubes. The IRAM-30m cubes 
are displayed in black, the interferometric cubes at 1, 3, and 6" are displayed 
in blue, green, and red, respectively. Top: Hybrid synthesis or IRAM-30m 
data cubes. Middle: PdBI-only data cubes or the 6'' extended data cube 
smoothed to 22.5^^ and subtracted from the IRAM-30m cube. Bottom: PdBI- 
only cubes subtracted from the hybrid synthesis cubes or the 6" PdBI-only 
data cube smoothed to 22.5^^ and subtracted from the IRAM-30m cube. Left: 
Histograms of absolute noise values. Right: Histograms of noise normalized 
by the median noise value. 
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Figure 12. Cumulative histograms of the signal above the 5cr-level for the 
^^CO (1-0) (dashed line) and ^^CO (1-0) (solid lines) cubes. The IRAM- 
30m cubes are displayed in black, the interferometric cubes at 1, 3, and ^" 
are displayed in blue, green, and red, respectively. Top: Hybrid synthesis or 
IRAM-30m data cubes. Middle: PdBI-only data cubes or the ^" extended 
data cube smoothed to 22.5^^ and subtracted from the IRAM-30m cube. Bot- 
tom: PdBI-only cubes subtracted from the hybrid synthesis cubes or the ^" 
PdBI-only data cube smoothed to 22.5^^ and subtracted from the IRAM-30m 
cube. Left: Histograms of raw brightnesses. Right: Histograms of bright- 
nesses normalized to the maximum brightness. 




The brightness noise levels of the hybrid synthesis data 
cubes are mostly set by the interferometric radiometric noise 
because the increase of integration time from the IRAM-30m 
to the PdBI does not match the gain in angular resolution (see 
Sect. |2.3.2| ). Hence, the median rms noise (see Table |4]) and 
the noise cumulative histograms (see Fig.[TT]) are almost iden- 
tical for the PdBI-only and hybrid synthesis data cubes. On 
the other hand, the brightness noise levels for the subtracted 
data cubes are mostly set by the single-dish data. Indeed, we 
subtracted two data sets whose noise is partially correlated: 
The noise coming from the PdBI visibilities is common. Sub- 
tracting both data cubes result in noise properties close to the 
IRAM-30m noise properties. In particular, the median rms 
noise in the subtracted data cube is more than one order of 
magnitude smaller than the median rms noise of the hybrid 
synthesis cube at V resolution. 

Figure [12] shows the cumulative histograms of the cube 
brightness using a similar layout as Fig.fU] In the right col- 
umn, the histograms are normalized by the maximum bright- 



ness. The histograms were computed using the full PAWS 
field of view but only the parts above the 5 a brightness noise 
level are displayed. The hybrid synthesis and PdBI-only his- 
tograms are similar both qualitatively and quantitatively. The 
histograms are displaced towards higher brightness values 
when the angular resolution increases, implying that all the 
structures above the 5 a noise levels experience large beam 
dilution effects. For example, the maximum brightness in- 
creases by more than one order of magnitude from 1.3 to 
16K, when increasing the resolution from 22.5 to V . On the 
other hand, the histograms of the subtracted cubes are identi- 
cal within the noise constraints (even for the normalized his- 
togram of the \" cube as the value of the maximum bright- 
ness is relatively uncertain). Indeed, the maximum brightness 
of the extended component evolves only from 0.97 ± 0.07 to 
1.36 ± 0.07 K from 22.5'' (the IRAM-30m resolution) to ^" 
(the PAWS resolution for which the extended emission is best 
defined). This implies that beam dilution is negligible for the 
subtracted emission. This is consistent with this emission be- 
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ing structured mostly at spatial scales larger than 36'^ 

The (/,m,v)- volume filling factors and the median bright- 
ness values were computed on all the cube pixels with de- 
tected signal, i.e. their brightness is larger than 5 times their 
noise. Using this definition, the brightness of the compact 
component (measured on the I" PdBI-only cube) ranges from 
2.0 to 16.0 K, with a median value of 2.5 K and it fills less than 
2% of the (/,m,v)- volume. The brightness of the extended 
component (measured on the &^ subtracted cube) ranges from 
0.07 to 1.36K, with a median value of 0.14K and it fills 
~ 30% of the (/,m,v)-volume. 

The comparison of the noise and signal distributions of the 
cubes at different angular resolutions makes it clear why the 
deconvolution of the PdBI-only data recovers more flux at 3^' 
than V^ but almost the same flux at &^ and 3'^ Indeed, the 
subtracted emission has a median brightness of 0.14 K, which 
is ~ 0.35 times the hybrid synthesis brightness noise level at 
I" but ~ 1 .3 and 4 times the hybrid synthesis noise level at 3" 
and 6", respectively. In other words, any emission present in 
the PdBI-only data is only recovered when its signal-to-noise 
is large enough (typically > 5). 

3.2.2. Spatial distribution of the line moments 

Figure [13] summarizes the properties of the decomposition 
of the totalemission into the compact and the extended emis- 
sion. The spatial distributions of the peak temperature, line 
integrated emission (the 0th order moment), centroid veloc- 
ity (the 1st order moment), line full width at half maximum 
(computed as 2.35 times the line 2nd order moment), and the 
noise are presented from top to bottom. The total and compact 
emissions are presented at the best possible PAWS resolution 
(V^), while the extended emission is displayed at the resolu 



tion where it is best measured, i.e., 6". Figures |32j [33| [34 
and [35] (available in the electronic version only) show tEe 
same^composition at fixed angular resolutions from 1, 3, 6, 
and 22. 5 '^ respectively. This demonstrates how the different 
moments of each component of the emission vary with spatial 
resolution. 

The deconvolved intensity distribution is corrected for pri- 
mary beam attenuation, which makes the noise level spatially 
inhomogeneous. In particular, the noise strongly increases 
near the edges of the field of view (see, e.g., panels [m] to 
[o] of Fig. [13]). To limit this effect, the deconvolved mosaic 
is truncateoat its edge, giving an almost parallelogram field 
of view of ~ 10.5 square arcminutes. A comparison with the 
IRAM-30m data (Fig. |35ja-f]) shows that signal is present 
near or slightly beyond tEe edges of the PAWS field of view. 
This implies that the signal at the edge of the PAWS field of 
view is probably less well deconvolved from the contribution 
of emission outside the PAWS field of view. This effect is 
more important at the northern edge and to a lesser extent at 
the southern edge than at the western and eastern edges. 

Two kinds of artifact appear in the peak temperature and 
integrated intensity maps of the I" subtracted cube, which 
shows the extended component. First, a moire effect due to 
the undersampling of the field pointings in the mosaic appears 
as a slight modulation of the intensity at a typical spatial scale 
of ~ 1 8^' (see Fig.[32tf]). This is due to power aliasing in the 
uv plane ( Pety & Rodriguez-Fernan dez 2010). Second, the 
chicken-pox aspect at a spatial scale close to the synthesized 
resolution is a known artifact of the deconvolution method 



ments at 3" and 6", where the impact of these artifacts be- 
comes negligible. 

The subtracted cube reveals extended emission whose peak 
temperature distribution is barely detected in the V^ hybrid 
synthesis cube, as the signal-to-noise ratio of this emission 
ranges between 0.2 and 3.5 (see Fig.[T3ja-c]). This is why a 
relatively complex dedicated masking technique was devised 
to compute meaningful 1st and 2nd order moments for the V^ 
hybrid synthesis cube (see Secti on [B] for a detailed descrip- 
tion). The peak temperature (Fig.|13[c]) and integrated emis- 
sion (Fig.[l3jf]) maps are maximum along the major axis of 
the galaxy. This is expected because emission in a given ve- 
locity channel extends over a large 2D area near the major 
axis, while is mostly extended in one spatial dimension along 
the minor axis. Hence, the interferometer will recover "ex- 
tended" emission along the minor axis much better than along 
the major axis. This projection effect thus minimizes emis- 
sion along the minor axis in the subtracted cube. However, 
there is more than a major axis trend in the subtracted cube. 
Indeed, Figure [30] and ^l\ (available in the electronic version 
only) show an overlay of the signal-to-noise contours of the 
V^ PdBI-only data cube onto the signal of the &^ subtracted 
cube for a set of channels at negative and positive velocities. 
These figures suggest that the extended emission fills the cen- 
tral 55'^ bounded by the inner edge of the spiral arms, and 
then falls on the convex side of the arms at larger radii (out to 
- 85^0. 

While the peak temperature map exhibits a symmetric spa- 
tial distribution relative to the galaxy center, the integrated 
emission peaks in the southern part. Extended emission is 
completely absent in the 1' x 1 '-areas located roughly south- 
west (-60'',-20'0 and northeast (+60'',+40'0 of the center. 
The maps of centroid velocity indicate differences between 
the kinematics of the compact and extended emission. This is 
best seen when following the Okms~^-isovelocity line, i.e., 
NGC5194's systemic velocity, on Fig. [34jh and i] or on 
Fig. [35th and i]. The linewidth of the extended component 
(Fig.|13[l]) is largest inside a central circle of ~ 35'^ radius. 
LinewlShs are on average much larger for the extended than 
for the compact emission. The clearest exception is at the 
galaxy center, i.e., at radii smaller than 2.5'^ where the com- 
pact emission has high peak temperature, large linewidth, and 
large integrated emission. This is reminiscent of the proper- 
ties of molecular gas in the inner 180pc of our Galax y ("Morris 
& Serabyn 1996). Alternatively, [Kohno et"aL] ( |1996| ) andMat- 
[sushita et al.| ( |2004| ) interpret these as gas being entrained by 
the AGN radio jet. 

We verified that the large linewidth of the extended compo- 
nent is not caused by the contributio n fro m the error beams. 
Details are provided in the Appendix |C. 3 [ 



3.2.3. Azimuthal averages 



that we used (see Fig. 32 c]). The overall spatial repartition of 
the extended component is nevertheless correct as evidenced 
by the comparison with the spatial distributions of the mo- 



Figure [14] shows the azimuthal average (and the associated 
dispersionjaround the kinematic center of the deprojected im- 
ages of the peak temperature, the integrated emission, the ro- 
tational velocity, the modulus of the centroid velocity gradi- 
ent, and the line full width at half maximum computed for the 
V^ PdBI-only (blue curves), the V^ subtracted (red curves), 
and the 22. 5^' IRAM-30m (green curves) cubes. 

For the compact emission (PdBI-only cube), the inner 5'^ 
clearly display high peak temperatures and large FWHMs, 
implying large integrated line emissions. The molecular ring 
dominates from ~ 10 to 40'^ where the integrated emission 
and the peak temperatures are larger than in the disk (radii 
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Figure 13. Comparison of the spatial distribution (from top to bottom) of the peak intensity, integrated intensity, centroid velocity, the line full width at half 
maximum (i.e., 2.35 times the standard deviation in velocity), and rms noise of the ^^CO (1-0) emission for the hybrid synthesis (PdBI + 30m, left column) and 
the PdBI-only (middle column) cubes at ~ I" resolution, and the 6^^ -resolution subtracted cube (right column). The angular resolution is indicated by a circle in 
the bottom left corner of each panel. The intensity scale is shown on the right-hand side of each panel. The major and minor axes are displayed as perpendicular 
dotted lines. The dotted circles show the two inner corotation resonances at radii equal to 23^^ and 55^^ while the dashed circle indicates the start of the material 
arms at a radius equal to 85^^ (Meidt et al.|2013jf . 



larger than ~ 40^0- The peak temperature seems to increase 
from the ring to the outer disk, while the integrated emission 
stays mostly constant outside r ~ 80'^ On the other hand, 
the velocity FWHM decreases slightly from ~ 25kms~^ at a 
radius of ~ 30^' to 20km s~^ at radii larger than ~ 100'^ 

The extended emission (subtracted cube) has a typical 
peak temperature of 0.75 K in the central region and 0.5 K 
in the disk. The FWHM decreases by a factor of 2 from 
^ lOOkms"^ at 0'' to ^ 45kms"^ at 50'' and it then varies 
between 40 and 50kms~^ in the outer disk. Both properties 
result in a regular decrease of the integrated emission from 
^ 45Kkms"^ at the center to -- 15Kkms"^ at 50''. It then 
varies between 10 and 15Kkms~^ 

The compact emission has on average a peak temperature 
twice as large as the extended emission. In contrast, the ex- 



tended emission has a velocity FWHM at least twice as large 
as the compact emission (except near the center). Both effects 
almost compensate to yield similar integrated line emissions 
for both components. 

The dispersion of the peak temperature and of the integrated 
emission is larger for the compact emission than for the ex- 
tended one. This is a consequence of the fact that the compact 
emission is structured at all scales down to or below the an- 
gular resolution while the filtered flux is structured mostly at 
scales larger than 36". The dispersion of the FWHM measure- 
ment is similar in both the compact and extended emissions. 
Indeed, its azimuthal average is computed only where there is 
enough signal to define it (i.e., on a small fraction of 360°), 
while the peak temperature and integrated emission are aver- 
aged over 360° (at least up to a radius of 85"). 
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Figure 14. Deprojected azimuthal averages around the kinematic center for 
the IRAM-30m (green curves), the PdBI-only (blue curves) and the subtrac- 
tion of the PdBI-only from the hybrid synthesis (red curves) cubes. The solid 
lines display the averages while the dashed lines give the azimuthal averages 
plus/minus the azimuthal standard deviations. From top to bottom, the panels 
present the peak temperature, line integrated emission, the rotational velocity, 
the modulus of the centroid velocity gradient, and 2.35 times the line 2nd or- 
der moment as a function of radius. The vertical dashed lines indicate the two 
inner corotation resonances at radii equal to 23^^ and 55^^ while the vertical 
dotted line show the start of the material arms at a radius equal to 85^^ (Meidt 
|et al. 201 3). Th e radial zone where the averages are affected by edge eftects 
(see Sect. |3.2.3) is highlighted in grey. 



The rotation velocity of each component was measured by 
fitting tihed rings with fixed systemic velocity to the line-of- 
sight velocity field using the GIPSY task ROTCUR. In both 
cases, we assume the kinematic center listed in Table [T] and 
we adopt a constant position angle (173°) and inclination 
(21°), as estimated from the more radially extended Hi emis- 
sion mapped at lower r esolution by the THINGS project (see 
[Colombo et al.||2013b| for more details). The middle black 
curve, labeled "Model", is a 3-parameter fit of the measured 
rotation curve. The inner part of this fit (inside lOO'O com- 
pares well with what we would expect if the stars (traced at 
3.6 /im) dominate the baryonic mass (Meidt et al. 2013). 

The rotational velocity of the extended emission (red curve) 
is increasing almost monotonically with radius, while it oscil- 
lates twice for the compact emission (blue curve) as an ef- 
fect of the streaming motions and corotation resonances (and 
not the bulge). Moreover, the centroid velocity of the ex- 
tended emission is typically closer to NGC5194's systemic 
velocity by 50kms~^ at radii smaller than 35^' where an inner 
stellar bar dominates the dynamics. At larger radii, it over- 
laps the rotational velocity curve of the compact emission. 
The modulus of the centroid velocity gradient is around 6 
and 10kms~V'' f^^ the extended and compact components, 
respectively. Hence, the kinematics of the extended emis- 
sion vary much more smoothly on the plane of the sky than 
the kinematics of the compact emission, which is very much 
affec ted by streaming motions (This is further discussed in 
Sect. [333] ). 

Any intrinsic behavior beyond a radius of 85 '^ must be inter- 
preted with caution as the azimuthal averages reach the edges 
of the field of view in its smallest dimension. Two effects hap- 
pen: 1) The noise increases sharply at the mosaic edges (see 
the bottom left panel of Fig. [13]), and 2) the azimuthal aver- 
ages miss the outside interarm regions, which occupy a larger 
and larger fraction of the area as the distance from the center 
increases. However, the comparison of the averages of the ex- 
tended and compact emission at each radius is meaningful as 
the averages are made on the same ellipse portions. 

3.2.4. Kinematics along the major axis 



Figure 15 compares the position- velocity diagrams along 
the major axis of the compact emission (PdBI-only, green 
contours) and the extended (grey image) emission. The red 
curve is the measured rotation curve. It matches well the over- 
all velocity variation along the major axis. The middle blue 
curve is the 3-parameter fit of the measured rotation curve, im- 
plying an overall inclination of the galaxy on the plane of sky 
of 21 ± 3° (see Sect. [3.2.3] ). The two other blue curves show 
the same velocity model for inclinations of 15° and 27°, in 
order to give an indication of the effect of inclination on the 
kinematics. 

This diagram confirms that the linewidth is much larger 
for the extended than for the compact emission (with the 
possible exception of the molecular gas at the center of the 
galaxy). The extended emission has a parallelogram shape 
with gas emitting at forbidden velocity in the [-35'',+35''] ra- 
dius ran ge. This is a typical s ignature o f nuclear bar kinemat- 
ics (g.g. ,[Binney et al.|[1991[[Garcia-Burillo & Guelin 1995; 



[Garcia- Burillo et al. 1999]k the distribution of the extended 
emission in PV-space might in addition or otherwise indicate 
that it lags the compact emission. The parallelogram shape 
is not symmetric, and emission is absent in the region near 
(+1 5^^,-25 km s~^). At and inside this position (and its mirror. 
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Figure 15. Position- velocity diagrams of the extended emission (hybrid synthesis minus PdBI-only) imaged at 6" along 3 axes at PA = 173° (major axis, central 
panel), and 173 ± 6° (left and right panels) shown in grey-scale. The scale is shown in the right-hand side (K [rmb]). The green contours (at levels 2, 4, 6, 8, 10 
K [T^hV) display the hybrid synthesis emission at \" resolution. The red curve is the azimuthally averaged rotational velocity profile and the central blue curve 
is a smooth 3-parameters fit of the red curve with a fixed inclination of 21°. The two other blue curves diplay the same model but for an inclination of 15° and 
27° . The vertical lines show the two inner corotation resonances at a radius of 23^^ and 55^^ ( JMeidt et al.|2013) . 

near -15^0. emission in the extended component that falls be- 
low the rotational velocities exhibited by the compact emis- 
sion (red curve) might arise with a genuine lag. 

3.3. Interpretation 1: Two CO disks - Thin and thick 

After an intermediate summary of the main observational 
properties of the two CO emission components, we will trans- 
late them into physical properties of the gas traced by the 
^^CO (1-0) emission. To do this, we will first summarize the 
expressions for the gas scale height, mid-plane pressure and 
gas density as a function of the gas and stellar surface den- 
sities and vertical velocity dispersions (Koyama & Ostriker 
|2009). We will then discuss how to apply these expressions 
to M5 1 . In particular, we will show how the contribution of 
the streaming motions to the CO linewidth can be estimated 
to derive an accurate vertical gas velocity dispersion. 
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3.3.1. Intermediate summary and consistency checks 

The emission filtered out by the PdBI interferometer ac- 
counts for about half of the total flux imaged in the hy- 
brid (PdBI+30m) synthesis. The subtraction of the PdBI- 
only from the hybrid synthesis cubes shows ^^CO (1-0) emis- 
sion mostly structured at angular scales larger than ?>6" , i.e., 
~ 1.3 kpc. Its brightness temperature ranges from 0.1 to 1.4 K 
with a median value of 0.14K. It covers about 30% of the 
PAWS field of view. Its spatial distribution surrounds the 
bright spiral arms. While the integrated line emission peaks 
at 48Kkms~^ near the major axis in the south-western galaxy 
quadrant, the peak brightness along the major axis ranges 
from 0.7 to 1.4K in the northeast and 0.5 to 0.8 K in the 
south-western quadrant. This emission is thus faint and ex- 
tended. In contrast, the emission observed by the interferom- 
eter is compact and bright. Indeed, its brightness temperature 
ranges from 2 to 16 K with a median value of 2.5 K and it 
covers less than 2% of the PAWS field of view. Rotational 
velocities estimated from the centroid velocity map of the ex- 
tended component are closer by 50- lOOkms"^ to the sys- 
temic velocity than the ones of the compact component inside 
a circle of 35'^ radius. Outside this radius, the rotational ve- 
locities of both components have the same order of magni- 
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Figure 16. Dual Gaussian decomposition of the THINGS HI spectrum aver- 
aged over the PAWS field of view after having aligned all the individual spec- 
tra along the velocity axis according to their centroid velocity value. Top: 
Brightness of the average spectra as a black histogram. The dual Gaussian 
fit is displayed as the green line and the two individual Gaussians are shown 
in red and blue. Bottom: Residual brightness after subtraction of the dual 
Gaussian fit. 





Table 6 

Results of the dual Gaussian fit (see Fig 


m- 


Gaussian Tpeak 

# (K[7;nb]) 


Velocity 
(kms-i) 


Width 
(kms-i) 


Area 
(Kkms-i) 


1 

2 


5.7 
7.2 


-1.4±0.4 

-f2.5±0.2 


80±2 
34±1 


490 ±13 
260 ±14 



tude. But, the rotational velocity curve of the compact com- 
ponent oscillates around the modeled velocities, while the ro- 
tational velocity curve of the extended component smoothly 
increases and mostly lies below the modeled velocities. The 
line FWHM of the extended component is twice as large as 
that of the compact component. 

We made two additional consistency checks. First, 
C PROPS decomposes the \" hybrid synthesis cube into two 
different components: 1) "Clouds" which account for 55% of 
the total flux, and 2) "Intercloud" gas which "surrounds" the 
GMCs and accounts for the remaining flux ( [Colombo et al. 
2013a). These numbers remain stable when the analysis is 
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done using the ^" hybrid synthesis cube. This decomposition 
result reflects the fact that we have a faint extended compo- 
nent in addition to the bright compact ^^CO emission. Second, 
we checked whether existing Hi observations are consistent 
with the possibiHty of having a narrow and broad Hnewidth 
component. We started from the THINGS Hi cu be imaged at 
11.9^^ X 10. C resolution with natural weighting (Walter et al. 
20081 ), as it maximize s the signal-to-noise r atio. We applied 
the shuffle method of |Schruba et al.| ( [2011| ), i.e., we shifted 
each spectrum of the Hi cube to a common velocity scale 
by removing the systematic velocity field structure as mea- 
sured from the centroid velocity. The spectra were then aver- 
aged over the PAWS field of view. Only a dual Gaussian can 
accurately fit the Hi spectrum. Figure [16] displays the shuf- 
fled, averaged spectrum, its Gaussian decomposition and the 
residuals. Table [6] presents the quantitative results of the dual 
Gaussian fit. The total flux is divided approximately equally 
in both components while the FWHM is approximately twice 
as large in one of the components. This is consistent with a 
separation into two emission components with very different 
linewidths. 

3.3.2. Expressions for the gas scale heights, mid-plane pressures, 
and gas densities 

Using the vertical momentum and Pois son equations aver- 
aged over the horizontal plane of the galaxy, [Koyama & Qs- 
triker| ( [200 9) obtained a second-order differential equation for 
the averaged vertical density profile. Solving it, they showed 
the averaged gas density (p) and pressure (P) are approxi- 
mately given by Gaussian profiles of the height z, i.e., 

p(z) = po exp ( -^ j and P(z) = Po exp ( -^ j • (1) 
In these equations, the gas scale height H is given by 
H= -^ 



V47rGp^A + VA2 + r 



(2) 



where a^ is the thermal plus turbulent velocity dispersion per- 
pendicular to the galactic disk, p^ is the stellar density, G is 
the gravitational constant, and A is a dimensionless factor that 
measures the relative densities of the gaseous and stellar disks. 
It can be expressed as 



A = 



2p^cr2' 



(3) 



where Sgas is the gas surface density. Once the gas vertical 
scale height is known, the gas mid-plane density and pressure 
can easily be derived with 



Po = 



27rH 



and Po = a,po. 



(4) 



These expressions for the gas scale height, mid-plane pres- 
sure, and mid-plane density take into account 1) gravity forces 
that both the stars and gas exert, and 2) turbulent and thermal 
hydrodynamic pressures. However, they still are lower limits 
as they neglect any contribution from the magnetic field. 

Finally, for an isothermal, self-gravitating stellar disk, the 
stellar surface density, the stellar volume density, the stellar 
vertical scale height, H^, and the stellar vertical velocity dis- 
persion, cr^, are linked via 



P* = 



2HJ 



and H^ = 



err 



ttG^. 



(5) 



The A factor can then be rewritten as 

4 2. 

A ~ , 



(6) 



where Q^ and 2gas are the Toomre's gravitational stability pa- 
rameters for the stellar and gas disks. 

3.3.3. Application to M51 



Figure [17] displays how the previous expressions are ap- 
plied to the case of M5 1 as a function of galactocentric radius 
for the total gas at 22. 5 '^ resolution (green curves), the ex- 
tended component at 6" resolution (red curves) and the com- 
pact component at V^ resolution (blue curves). Radial zones 
where the results should be interpreted with caution are high- 
lighted in grey: 1) at radii larger than 85'^ the azimuthal av- 
erages start to reach the edges of the observed field of view; 
and 2) from to ~ 45'^ the assumption that the stars are dis- 
tributed in a disk breaks down, so that the stellar scale height 
and volume density are not well constrained (see below). 

The gas mass surface densities for each component are 
computed from the azimuthal averages of the ^^CO (1-0) in- 
tegrated emission using the Galactic value of the Xco factor 
and taking into account the presence of helium. The stellar 
surface density is derived from the 3.6 /im emission ( |Meidt 
let al. 2012). The stellar velocity is computed following Bot 
[tema ( 1993 ) and Boissier et al. (2003 ) who showed that, for a' 
flat exponential disk, it falls off from the central dispersion ctq 
according to 

/ r \ 

(7) 



a^ = ctq exp - 



2He 

where r is the galaxy radius and Hb is the disk scale length of 
the 5-band. jMcElroy (1995) estimates the central stellar ve 

"^ and Trewhella 



locity dispe rsion in M51 to be ctq = 113kms 



etal 



estimates that Hb = 2.82kpc. This gives us an es- 

timate for tfie stellar scale height and volume density, accord- 
ing to Eq. [5] We emphasize that this assumes an isothermal, 
self-gravitatmg stellar disk. This assumption clearly breaks 
down in the center of M5 1 , where the bulge and nuclear bar 
dominate, i.e., at radii less than ^ I'S" where the 3.6 /im sur- 
face brightness profile exhibits steepening, and near the lo- 
cation of the ba r corotation radius as estimated from gravi- 
tational torques (jMeidt et al.||20131). In fact, the application 
of Eq. [5] results in a stellar scale height that decreases with 
radius irom ^ \'S" to the galaxy center. We therefore opted 
to adopt a constant, lower limit to the scale height inside this 
zone by extrapolating the value given by Eq. |5] at a radius of 
45 '^ inward. 

The vertical velocity dispersion of the gas can be estimated 
as the line 2nd order moment for a face-on galaxy. The in- 
clination of M51 onto the line of sight is small but non-zero, 
implying that the line 2nd order moment is only a first order 
approximation. Indeed, the systematic motions averaged in- 
side the beam of the observations contribute to the line 2nd 
order moment. This is more problematic in the case of M5 1 
because the streaming motions are known to be large for this 
galaxy. Appendix O shows that the vertical velocity disper- 
sion can be estimated as 



^z ^ ((^obs-Vcent) )" 



n2 



Igrad(Vcent) 



2.35 



(8) 



where () symbolizes the brightness- weighted average over the 
line profile, Vobs is the velocity projected along the line of 
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Figure 17. Azimuthal averages around the kinematic center, from top to bottom: Mass surface densities, vertical velocity standard deviations {i.e., line full 
width at half maximum divided by 2.35), gas-to-stellar density ratios (A factor), gas thermal -i- turbulent mid-plane pressures, gas mid-plane densities, vertical 
scale heights for the stars (pink curves), the IRAM-30m (green curves), the PdBI-only (blue curves) and the subtracted (red curves) cubes. The bottom panel 
displays the ratios of the scale heights (black) and mid-plane densities (cyan) computed from the PdBI-only and subtracted data sets. The horizontal lines indicates 
typical values of the different parameters for the compact and extended components. The vertical dashed lines indicate the two inner corotation resonances at 
radii equal to 23" and 55^^ while the vertical dotted li ne show the start of the material arms at a radius equal to 85^^ |Meidt et al.|2013| . The radial zones where 
the results should be interpreted with caution (see Sect. |3.3.3) are highlighted in grey. 
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PdBI-only 



Hybrid synthesis - PdBI-only 




20 40 

|grad(v,,,t)l-(0/2.35) [km.s"'] 



Figure 18. Joint distributions of the line 2nd order moment as a function of 
the centroid velocity gradient modulus times the half primary beam width for 
the 30m (left), PdBI-only (middle) and subtracted (right) cubes. The lines 
have a slope of 0.2 (green, soHd), 0.5 (green, dashed), 1 (black, soHd), 2 (red, 
dashed), and 5 (red, solid). 



sight, Vcent the Hne centroid velocity, |grad(Vcent)| the modu- 
lus of the gradient of the centroid velocity, and 6 the resolu- 
tion beamwidth of the observations. In this equation, the first 
term is the square of the second moment and the second term 
estimates the contribution of unresolved systematic motions. 

Figure [18] diplays the joint distributions of these two quan- 
tities for tEe IRAM-30m cube and the l" PdBI-only and sub- 
tracted cubes. At the resolution of the IRAM-30m data, the 
unresolved systematic motions contribute significantly to the 
value of the second moment. This behavior is clearly split 
in the decomposition between compact and extended com- 
ponents at \" . The unresolved systematic motions are neg- 
ligible for the extended component. For the compact com- 
ponent, they contribute to less than 34% of the linewidth for 
50% of the data. Hence, the streaming motions are seen in 
the compact component but they are not seen in the extended 
component. We stress that the subtracted cube measures the 
extended component at an angular resolution of \" because 
it results from the subtraction of two cubes whose resolution 
is \" , namely the hybrid synthesis and the PdBI-only cubes. 
Beam smearing of unresolved systematic motions must there- 
fore be considered only at angular scales lower than I" . As a 
corollary, if the large linewidths are due to beam smearing of 
unresolved streaming motions, then the linewidths should de- 
crease when increasing the imaging angular resolution. This 
effect is observed in both the hybrid synthesis and PdBI-only 
cubes (panels [j] and [k] of Fig. [32] to [34]), while the sec- 
ond moment of the subtracted cube (panels [1] of the same 
figures) stays basically constant when increasing the angular 
resolution from 22. 5'^ to \" . We thus deduced that the large 
linewidths of the extended component are not caused by un- 
resolved streaming motions. 

As the unresolved systematic motion can be larger than the 
second moment for the compact component, we only used 
the second moments to compute the vertical velocity disper- 
sion, implying that the scale heights and mid-plane pressures 
derived are slightly more robust for the extended component 
than for the compact component and the IRAM-30m data. 

3.3.4. Results 

Using these inputs, we computed the stellar vertical scale 
height, and volume density. We also computed the gas-to- 
stellar density ratio (A factor), the gas mid-plane pressures, 
densities, and scale heights for all the molecular gas (using the 
IRAM-30m data), the compact component (using the PdBI- 
only data), and the extended component (using the subtracted 
cube). We neglect the contribution from atomic gas traced by 
Hi emission, as this gas represents only between 2.5 and 30% 
of the molecular gas mass for radii from to 130^' (see Fig. 42 



of |Leroy et al.|2008t [Schuster et al.|20Q7 l). 

Our estimates agree very well with expectations. For ex- 
ample, the A factor is less than 1 for the extended compo- 
nent, while it is equal to ~ 2.5 for the compact component 
for radii larger than 0.5 kpc. Using Eq. [6] this implies that 
Qcompact ^ 0-32extended, where Qcompact and ^extended are the 
Toomre factors for the compact and extended components, 
respectively. Said otherwise, the compact component is more 
likely to form stars than the extended component. The gas 
mid-plane pressure is about the same for both components and 
for the total gas. This probably reflects pressure equilibrium. 
From the mid-plane pressure, we estimate a molecular frac- 
tion Su /Efjj close to what is observed: Using the empirical 
formula which relates the molecular fraction Sfj /^Hl ^^ ^^^ 
mid-plane pressure ( [Blitz & Rosolowsky|2006| ) 



%9/%r 



4.3xl04Kcm-3 



0.92 



(9) 



we find a molecular fraction of ~ 26 and ~ 4 at radii equal 
to 1.0 an d 2.5 kpc, resp ectively, in good agreement with the 
results of |Leroy et al.| p008). Based on the high molecu- 
lar fraction, we thus expect a tight correlation between the 
star formation rate and the molecular gas surface density, but 
probably little relation to the a tomic gas surface density, as 
found by Kennicutt et al. ( 2007| ). 

The scale height of the extended component is almost con- 
stant, its value varying slighty betwen 190 and 250 pc for radii 
larger than 0.5 kpc. The scale height of the compact compo- 
nent varies more. It is ~ 40pc between ~ 1.5 and 3.5 kpc. 
It varies between ~ 20 and 40 pc beyond ~ 3.5 kpc and it 
decreases to ~ lOpc at r^ 0.5 kpc. In other words, the scale 
height of the compact emission is 5-6 times smaller than the 
scale height of the extended component from ~ 1 .5 to 3 .5 kpc. 
The compact component is between 5 and 20 times (typically 
10 times) denser than the extended component. At a radius of 
2.5 kpc, the extended and compact components have an aver- 
age volume density of ~ 1 and 10 H2 cm~^, respectively. For 
reference, the distribution of the molecular gas density in the 
Solar Neighborhood is ( ,Ferriere,2001 , ,Cox,2005^ 



n 



H? cm 



= 0.29exp[-(z/81pc)2], 



(10) 



i.e., a scale height of 57pc(= 81pc/>/2) and a mid-plane gas 
density of 0.29 H2 cm~^. The volume density is several orders 
of magnitude smaller than the expected densities of molecular 
gas. This is due to the fact that the molecular gas fills a small 
fraction of the galactic volume. 

3.4. Interpretation 2: A mixture of dense and diffuse gas 

The averaged volume density of the compact component is 
typically one order of magnitude larger than the one of the ex- 
tended component, pointing toward different kinds of molec- 
ular gas. In this section, we recall that 1) bright ^^CO (1-0) 
emission traces diffuse as well as dense gas, and 2) the value 
of the r(^^CO)/r(^^CO) ratio may be used to discriminate 
between dense and diffuse gas. We will then check this ratio 
forM51. 

3.4.1. Bright ^^CO (1-0) emission also traces diffuse gas 

Bright ^^CO (1-0) emission is generally associated with 
dense cold (typically n - 10^- lO^cm"^ and T - 10-20 K) 
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Figure 19. Joint distributions of the ^^CO (1-0) emission as a function of the ^^CO (1-0) emission. Top row: Brightness integrated over the line profile. Bottom 
row: Brightness in 5kms~^ channels. Contour levels are set to 2, 4, 8, ... 2048 and 8, 16, 32, ... 2048 points per pixel respectively for the top and bottom rows. 
First column: Full field-of-view (black contours). Second column: Radii below 35^^ (red contours). Third column: Radii between 35^^ and 150^^ (green 
contours). Fourth column: Radii larger than 150^^ (blue contours). The 3 same straight lines display 3 different ^^CO/^^CO emission ratios on each panel: 7 
(red), 8.3 (green) and, 11 (blue). 



molecular gas, where all hydrogen is molecular and all car- 



bon is locke d in CO. However, |Pety et al.| ( |2008| ) and |Liszt 
et al.| ( |2009t found surprisingly bright ^^ CO (1-0) lines (up to 
~ lOK) m the nearby environment of Galactic diffuse lines 
of sight (Ay '^1), where the hydrogen is partly atomic and 
partly molecular, and where the carb on occurs mostly in ion 
izedform, i.e., C^. 



higher density gas. In diffuse gas, ^iszt & Lucas ( 2000 ) and 
Liszt] (|2007) showed that if the amount of HCO"^ observed in 



the diffuse gas is fixed as a model parameter, it is easy to get 
large amount of CO in UV illuminated gas through electron 
recombination of HCO"^ 



Liszt et al. 



( 2010t explain such large ^^CO 
brightnesses in diffuse warm gas (typically n ~ 100-500cm~^ 
and r ~ 50- lOOK) by the fact that the gas is subthermally 
excited gas. I ndeed, large velocity gradient radiative trans- 
fer ni ethods ( Goldreich & Kwan||1974t [Scoville & Solomon] 
1974| show that 1) Wco/Nqo is large because of weak CO 
excitation in warm gas (50-100 K), and 2) Wqo oc Nqo until 
the opacity is so large that the transition approaches thermal- 
ization. Hence, relatively bright ^^CO lines may at least trace 
either diffuse warm or dense cold molecular gas. 

This is surprising because it is often argued that CO can- 
not survive outside dense molecular gas since chemical mod- 
els predict that several magnitude of visual extinction are 
required so that CO survives photo-dissociation. However, 
many absorption measures in the UV and millimeter domains 
show that CO is present in gas w hose hyd rogen column den 

sity is as low as 

Sonnentrucker et aLl2007| ). 

CO ch emistry in"5iifuse ga s is still far from being under- 
stood. [Sonnentrucker et aL] (|2007) found that a plot of the 
log(Ato) as a function of the log(A^fj ) can only be cor- 
rectly fitted by two power-law relationships, with a break at 
(A^co = 1.3 X 10l4cm-^%2 = 2.5 x lO^^cm"^), correspond- 
ing to a change in the production route for CO. The produc- 
tion routes of CO are well understood only in the regime of 



HCO+ + e"^CO + H. 



(11) 



10^2 cm-2 dLiszt 2008; S heffer et al ||2008 
TKe key point here is that the 



[Visser et al.| ( [2009| ) later confirmed this result. The next (still 
unsolved) question is how large quantities of HCO"^ form in 
the diffuse gas. 

3.4.2. The value of the T i}^ CO) / T i}^ CO) ratio may discriminate 
between dense and diffuse gas 

As bright CO (1-0) lines may be associated with diffuse or 
dense gas, this line alone ca nnot be used to differentiate both 
scenarios. [Lisztet al.|(2010t argue that the TC^C0)/TC^C0) 
ratio can discriminate dittuse and dense gas. Indeed, from 
wide (-a5^V^C0and ^^CO (1-0) maps of the Galactic 
plane, Polk et al1(|1988| ) measured an average ratio of line in- 
tensities R = TC^CO)/tC^CO) of 6.7 ± 0.7. Matching the 
beam areas at the two frequencies introduces an upward cor- 
rection factor of 1.1, /.^., 7? = 7.5 ±0.8. This value is higher 
than the typical factor measured f or the core of Giant M olec- 
ular Clouds, i.e., ~ 3-5 (see e.g. Frerking et al.||1982| ). Us- 
ing the observational fact that r(^ -^CO)/r(^^^) ^ 10-20 
in diffuse molecular gas (see e.g. |Knapp & Bowers||l988| ), 
|Polk et aL] ( [1988| ) deduced that this diffuse gas component 
significantly contributes to the large-scale ^^CO em ission of 
the Galaxy. More recently, [Goldsmith et al.| ( |2008| ) deduced 
from ^^CO and ^^CO wide-field mapping of the Taurus giant 
molecular cloud that about half the mass of the gas traced by 
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the ^^CO (1-0) emission comes from diffuse gas. 

The radiative and chemical properties of CO explain the 
large values of this ratio in diffuse molecular gas. From mea- 
surements of CO absorption lines against extra-galactic con- 
tinuum background sources, |Liszt & Lucas | (1998) showed 
that in local diffuse clouds, the ratio of ^"^CO and ^^CO col- 
umn densities is in the range 15 < NC^C0)/NC^C0) < 54, 
with this ratio declining with increasing N^^CO). Hence, this 
ratio usually differs strongly from the local interstellar ratio 
of the C elemental abundances, typically ^^C/^^C = 60. The 
value of the A^(^^CO)/A^(^^CO) ratio is due to the competi- 
tion between the fractionation of CO with C"^ and selectiv e 
photo-dissociation of the two CO isotopologues ( |Liszt|200"7] ). 



al., in prep.). The comparison of this map with the 22. 5 '^ 
decomposition between compact and extended sources (see 
Fig. [35]) will probably shed light on the origin of the [Cll]. 



Fractionation of CO with C"^, i.e.. 



i^CO + ^^C+^^^CO + ^^C^, 



(12) 



enriches CO in ^^CO. But selective photodissociation (i.e., the 
fact that ^^CO better self-shields from UV illumination than 
^^CO) more than counters this effect in diffuse gas. 

3.4.3. InM51 

Figure |3] indicates that the morphology of the ^^CO and 
^^CO moments are strikingly similar where both tracers are 
detected. To quantify this. Fig. [19] shows bi-dimensional his- 
tograms of the integrated intensity (top row) and brightness 
(bottom row) of the ^^CO (1-0) versus ^^CO (1-0) emission. 
Both sets of histograms display similar linear relationships 
whose width is typically related to the noise levels of both 
tracers. The left column displays the histogram computed for 
the full field of view, while the three other columns show how 
this histogram decomposes as a function of radius. For ref- 
erence, we drew the same three lines on all histograms. The 
slope of these lines were chosen to follow the ridges of the 
brightness histograms in the different radial ranges. 

While the typical value for the ^^CO/^^CO ratios is about 
8, there is a slight but significant increase of this ratio from 
the inside to the outside of the galaxy: This ratio is typically 
7 (blue lines and histograms) for radii below 35'^ 8.3 (green 
lines and histograms) for radii between 35 and 150^' and 11 
(red lines and histograms) for radii beyond 150'^ For com- 
parison, this ratio increases from 4.6 to 10.0 across the Milky- 
way molecular ring, i.e. from O.SRq to Rq (LisztetaL1984). 
This was an early indication that molecular gas near the Solar 
Circle has a high proportion of diffuse material. 

At the angular resolution of the IRAM-30m (22.5'0, the 
r(^^CO)/r(^^CO) increases from 6 to 11 with a typical value 
of ~ 8 in M51. Assuming a Galactic [^^C]/[^^C] elemental 
ratio, the comparison with our Galaxy points towards the in- 
terpretation that a significant fraction of the CO emitting gas 
is diffuse. Is this visible in the emission of C"^? iNikola et al.l 
([2001) published the first map of the 158 /im [Cll] emission 
at 55'' resolution for M51 obtained with the Kuiper Airborne 
Observatory. The interpretation of this emission is difficult 
because it originates mostly from the warm ionized medium, 
though the cold neutral medium contributes significantly to 
the [Cll] emission. In the cold neutral medium, the density so- 
lution is degenerate, i.e., the medium could have either a low 
(- 100-300cm-^) or a high (10^ - lO^cm"^) density. They 
conclude that 'A large fraction of the overall [Cll] emission 
in M51 can originate in an underlying extended medium." 
A 12^' -resolution [Cii] maps was observed as part of a Her- 
schel/PACS guaranted time project (PI: C. Wilson, Parkin et 



3.5. Discussion 

Here we discuss the structure of the gas that emits the ex- 
tended component — in particular the fact that it could be 
substructured like drops in fog — and the amount of gas that 
is extra-planar. We then summarize results about extra-planar 
gas traced in Hi and CO. We finish with a discussion of the 
possible origin of extra-planar CO emission. 

3.5.1. Structure of the extended component 

Giant molecular clouds are often thought to be composed of 
small clumps that are unresolved, i.e. a set of point sources. 
Hence, the interferometer should recover all the flux of this 
set of point sources and short-spacings should not be needed 
for extra-galactic observations. This argument is incorrect for 
two reasons. The first reason is that the limited sensitivity 
of the interferometers limits the power of deconvolution to 
recover th e flux of point sources at low signal-to-noise ratio 
(See Sect. [2X2l and |2331 ). 

The second reason why the PdBI-only data may not recover 
the full flux of the source is more fundamental. We first as- 
sume that the source is a set of unresolved components filling 
a volume that projects onto a plane-of-sky area larger than 
the interferometer primary beam. If the unresolved compo- 
nents are typically separated by an angular distance larger 
than the synthesized beamwidth, the interferometer will in- 
deed recover all the flux because the GMC is observed as a 
set of separated point sources. In contrast, the interferometer 
will filter out most of the source flux when the typical angu- 
lar separation between the unresolved components is smaller 
than the synthesized beamwidth because the source emission 
appears like a flat source. Well known examples of this ef- 
fect are 1) fog which appears flat while made of water drops, 
and 2) a XlXth century pointilliste painting which appears flat 
when unresolved by the eye. For this reason, it is not possible 
to know before the observation of large complex sources like 
nearby galaxies whether the short-spacings will be needed, 
and multi-resolution CLEAN algorithms cannot help to solve 
this ambiguity. 

To explore this further, we speculate that the CO extended 
component is made of diffuse gas as found in the envelopes 
of giant molecular clouds. This gas would have a typical 
mid-plane density < lOOOcm"^, and a temperature < 200 K. 
In such conditi ons, th e ^^CO (1-0) is subthermally excited. 
Moreover, Sect. 1 3. 3. 4| indicates that the average volume den- 
sity is lcm~^, implying a typical volume filling factor of 
0.1%. The structure of this gas is probably filamentary, as 
observed in our Galaxy. 

3.5.2. Vertical mass distribution and extra-planar gas 

Figure |20| displays the vertical distribution of the molecular 
gas volume density, its decomposition into two components 
of different scale heights, and the percentage of mass above a 
given galactic height |z|. In order to sample a large fraction of 
the sensible parameter space, we computed 4 different cases: 
Two fractions of flux in the extended component (either 30 
or 50%) and two scale height ratios (5 and 10) between the 
compact and extended components. We see that the fraction 
of flux in both components has a significant impact on the 
gas distribution close to the galatic mid-plane. However, the 
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Figure 20. First and third columns: Vertical distributions of the total volume density (green), decomposed into the sum of two Gaussian of different scale 
height. The broad and narrow Gaussian are respectively called diffuse (red) and dense (blue). Second and fourth columns: Percentage of the vertical mass above 
a given altitude in scale height unit of the dense Gaussian for the total distributions plotted in the first and third columns. The vertical lines indicate the altitudes 
above which 1% (red), 3% (green), 10% (blue), 33% (cyan), and 100% (yellow) of the mass is located. The ratio of the integrated masses of the diffuse over the 
total (diffuse+dense) are 0.3 for the top row and 0.5 for the bottom row. The two right (respectively left) columns are for a diffuse scale height 5 (respectively 10) 
times higher than the dense scale height. This allows us to quantify the mass of molecular gas which is extra-planar in different scenarii. 



percentage of the total mass above a given height z is only 
marginally affected. In contrast, doubling the ratio of the scale 
height doubles the galaxy height above which a given mass is 
located. For instance, if we assume that the smallest scale 
height is 40 pc and a scale height ratio of 5, only 2% of the 
total mass is above 400 pc while this proportion increases to 
20% when the scale height ratio is 10. We thus estimate that 
between 2 and 20% of the molecular gas is extra-planar, i.e. it 
lies at a galactic height at least 10 times larger than the scale 
height of the dense molecular gas. 

3.5.3. Hlth ick disk 
Lagging, thick Hi layers have long been detected in external 



galaxies (see, e.g., |B oomsma et al. 2005 ; B arb ieri et al . |20Q5 
IQosterloo et~aLl|2QQ7t iBoomsma et al.| |2008 , in NGC 253, 
NGC 4559, NGC 891 and NGC 6946). The Hi gas in these 
galaxy halos amounts to 3-30% of the total Hi mass. For edge- 
on galaxies, the Hi emission extends up to 12-22 kpc from the 
g alactic mid-planes. 

[Miller et al.| ( |2009| ) studied the prototypical face-on spi- 
ral galaxy M83. Tney found a spatially extended compo- 
nent rotating in the same sense but 40-50 km s~^ more slowly 
in projection, with a line-of-sight velocity dispersion of 10- 
15kms~^ The spatially extended structures are coincident 
with the optical spiral arms. They interpreted this compo- 
nent as a vertically extended disk rotating in the same sense 
but about 100 km s~^ more slowly than the kinematically cold, 
thin disk. It contains 5.5% of the total Hi mass within the 
stellar disk of the galaxy. 

3.5.4. A CO thick disk in the edge-on NGC 891 galaxy 

Among the previous galaxies, NGC 891 is particularly in- 
tere sting. Its characteristics are very similar to the Milky Way 



and lGarcia-Burillo et al.| ( |1992| ) found in this system the first 
evidence for extra-planar CO emission in an edge-on galaxy. 
They detected molecular emission 1-1. 4 kpc above the disk. 



wide-spread along the major axis. The associated Gaussian 
scale height is > 600 pc, i.e. typically 2 to 3 times larger than 
that deduced here for M5 1 . This CO scale height is confirmed 
by a new IRAM/HERA map, which now images most of the 
thick disk of the galaxy (Garcia-Burillo et al., in prep.), and 
it is consis tent with the s cale height derived from the PAH 
emissi on (|Rand et al.|2008p. 

Garcia-Burillo et al.|(|1992|) estimate this halo emission rep- 
resents less than 20% of the molecular mass in the disk. How- 
ever, they also state that 40 to 60% of the thin disk emission 
must be in a low density component (diffuse gas) to explain 
the typical value (8.5) of the ratio of the ^^CO (1-0) to ^^CO 
(1-0) brightnesses. Hence, the halo molecular emission could 
be the tip of an independent component mixed with the dense 
molecular thin disk. In this case, the extended component 
would amount for a much larger fracti on of the entire molec- 
ular gas. Indeed, in a follow-up study Sofu e & Nakai| ( p^93] ) 
mapped the molecular gas at a distance of 3.5 kpc from NGC 
891 's galaxy center, i.e. , at the ring northern side. Using a dual 
Gaussian fit through the emission profile extracted perpendic- 
ular to the galaxy plane, they estimate 1) that the extended 
component has a typical scale height of 2.1 kpc and 2) that it 
could account for 50% of the entire molecular gas. 

3.5.5. A diffuse CO thick disk in our Galaxy? 

Mapping 3 strips at constant longitude within a latitude 
range of [-3.5°, +3.5°], Dame & Tha ddeus] fT994\ detected 
a thick molecular disk in the inner Galaxy about 3 times as 
wide as that of the dense central CO layer, and comparable in 
width to the thin Hi layer, i.e. the height profile can be fitted 
with a Gaussian FWHM of '^ 23 Ope. This data also suggests 
that the high latitude gas lies mainly above the spiral arms. 
The mass of this gas would be 15% of the total if it belongs to 
a distinct molecular component. 

Combining £'b-v reddening. Hi absorption and ^^CO emis- 
sion measur ements along many diffuse lines of sight in the 
Milky Way, [Liszt et al.| ( ,2010) recently found that the CO 
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luminosity per H2 molecule for diffuse gas is standard, i.e. 
Xco = 2.0 X 1 0^Qcm~^/(Kkms~ ^). This standard value was 
introduced by [Polk et al.| ( |1988] ) to take into acco unt the con- 
tribution of the diffuse gas to t he Xco factor. | Liszt et aL] 
12010) and [Liszt & Pe ty (2012) then deduced that the dif- 
fuse gas contribution to the total CO luminosity seen looking 
down on the Milky Way is 0.47Kkms~^ to be compared to 
0.75Kkms~\ the result from Galactic surveys. Hence, the 
contribution of diffuse gas to the total CO luminosity in our 
Galaxy would be between 39% (= 100 x 0.47/(0.47 + 0.75)) 
and 63% (= 100 x 0.47/0.75), depending on the fraction of the 
low brightness diffuse component that was already detected 
by the Galactic surveys. 

The typical sensitivity of wide-field CO surveys is 
IKkms"^ In addition, they have been mostly observed be- 
fore the advent of wide bandwidth receivers/spectrometers, 
which were installed in current facilities from ~ 2005 on- 
wards. The total velocity bandwidth of the well-known wide- 
field Galactic CO surveys is typically only ~ 300kms~^. Fi- 
nally, the extension of the diffuse component we detect in 
M51 is large, several hundred of parsecs. Hence, such a com- 
ponent in our Galaxy could have been confused with low level 
baselines in wide-field CO surveys. 

3.5.6. Possible origin of the extra-planar CO emission 

The following three scenari os are typic ally evoked to ex- 
plain gas outside the disk (e.g. [Putman et al. 2012): expelled 
gas from the disk via galactic fountains or chimneys, accre- 
tion from the inter-galactic medium and/or tidal debris from 
galaxy interactions. Except for the first scenario, it is expected 
that the involved component reaches very large (more than a 
few kpc) scale heights. 

Therefore we focus in the following discussion on observa- 
tional and theoretical evidence for a potential flow of dense 
cold material from the disk to the halo caused by massive star 
formation. The most extreme case to be considered might be 
molecular outflows driven by stellar winds or nuclear AGN, 
as seen, e.g., in CO in the nearby starburst galaxy M82 (^Wal-| 
ter et al. [2002|), and recentl y by Herschel in several nearby 
(U)LIRGs lSturm et al.f 201 1 ). However, typically these phe- 
nomena are restricted to the centers of galaxies. 

Evidence for dense material outside the thin disk comes 
from observations of distinct optical extincti on features in 
edge-on galaxies that trace extra-planar dust. |Howk| ( |2005| ) 
summarizes the findings on extra-planar dust in nearby galax- 
ies: extra-planar dust is found out to scale heights of z < 2 kpc, 
correlates well with regions of massive star formation in the 
disk and the total amount of star formation, and appears to 
reside in cold dense gas with densities of ^(H) > 25 cm~^ . In 
this context the detection of abundant extra-planar PAH emis- 
sion is interesting. Rand et al. (2011) find PAH emission out 
to ~ 0.5 - Ikpc above the disk plane in two edge-on spiral 
galaxies with detected extra-planar dust and detect emission 
from the mid-IR 17 /im H2 line out to distances of 2 kpc in 
their targets. They speculate that massive star formation in the 
disk is the cause for the extra-planar cold interstellar medium 
detected. 

The resolved CO emission in M51 is preferentially found 
along the convex side of the spiral arms where massive star 
formation is ongoing and has an inferred scale height of 
z ~ 200 pc. These two findings are very similar to results 
of the studies of extra-planar dust and PAH emission. Thus 
we speculate that indeed star formation via galactic fountains 



or chimneys has transported some of the molecular material 
away from the disk. 

Several simulations of disk galaxies with star formation 
also suggest that this explanation might be valid. The simula- 
tions have shown that cold (T < 200 K), dense gas can reach 
heights of 100-200pc above the plane, but only when stellar 
feedback is in cluded (Wada 2008; Koyama & Ostriker 20091 
Dobbs et al.1[2QllnAcreman et al.||2012^ |HiU et al.||20n] ). 
Dobbs et al.| ( |201l] ) find the scale height of the cold gas is 
^ 50- 1 00 pc, depending on the level of f eedback in the simu - 
lation. From the top panel of Figure 9 of Dobbs et al.| ( [2QTT| ), 
it is seen that gas at heights of a few 100 parsecs is around 
10-25 - IQ-^'^gcm-^ or 20-200H2 cm'^). Gas at these den- 
sities is not typically molecular. However, the top of Figure 
14 of Dobbs_etalj 02OO8 ) shows that gas which has already 
been dense and molecular can retain a high molecular frac- 
tion down to densities of ~ 2 H2 cm"^ before the H2 is pho- 
todissociated. Thus the simulations indicate that a possible 
explanation of an extended diffuse component is that stellar 
feedback pushes gas out to large distances above the plane, 
but this gas remains molecular. The likelihood of this de- 
pends on the gravitational contribution from the stars/gas in 
the vertical direction, the local chemistry of H2 formation and 
destruction, the effects of the stellar feedback and the surface 
density of the gas. 

4. SUMMARY AND CONCLUSIONS 

We described in detail the calibration and construction of 
the PAWS i^CO (1-0) imaging of the central - 10 x 6kpc 
(-- 270'' X 170'0 in M51, using observations from both the 
PdBI and 30m telescope. The achieved spatial resolution of 
40pc (l.l'O is close to the typical size of galactic GMCs, 
and at least 10 times smaller in area than previous in terfer- 
ometric maps obtained at OVRO, BIMA and CARMA ( |Aalto 
[eraL][T999l [Heifer et al.| [20031 [Koda et al.[[2QTT] ). THe 
median brightness sensitivity of 0.4 K in 5kms~^ channel 
spacing corresponds to 8.6 M© pc~^. The total flux in the 
PAWS field of view and in a velocity range of ±120kms~^ 
around the LSR systemic velocity is 64% of the total M5 1 
flux of 1.4 X 10^Kkms~^pc^, i.e., a molecular gas mass of 
6.2 X 10^ Mq (helium included). The mean CO integrated in- 
tensity and molecular mass surface density inside the PAWS 
field of view are 18Kkms~^ and 77M0pc~^, respectively. 

The interferometer recovers only (50 ± 10)% of the total 
flux. The remaining flux is mostly distributed on scales larger 
than 1.3 kpc (36'0. Hence, the flux is about equally distributed 
into a spatially extended and a spatially compact component. 
Using the hybrid synthesis (PdBI+30m) and the PdBI-only 
deconvolved results, we established that the extended compo- 
nent has the following properties: 1) It has a median bright- 
ness temperature of 0.14K, about 18 times fainter than the 
compact component. 2) It covers about 30% of the PAWS 
field of view, about 15 times as much as the compact compo- 
nent. 3) Its plane-of-sky kinematics evolves smoothly, while 
plane-of-sky kinematics of the compact component are more 
affected by streaming motions. 4) Its linewidth is typically 
twice as large as that of the compact component. 5) Outside 
the region dominated by the galaxy bulge and inner nuclear 
bar, its typical scale height is ~ 200 pc, or five times the com- 
pact component scale height. 6) Its typical gas mid-plane 
pressure is '^ (2-4) x lO^Kcm"^, in approximate pressure 
equilibrium with the compact component. 7) Its typical av- 
erage mid-plane gas volume density is ~ IH2 cm~^, 10 times 
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less dense than the compact component. 

We estimated between 2 and 20% of the total molecular 
mass to be extra-planar, i.e. at galactic heights larger than 
400 pc. We emphasized that, while the emission of the ex- 
tended component is mostly structured at spatial scale larger 
than 36'^ it is probably made of unresolved filamentary struc- 
tures, which could typically fill 0.1% of the volume of the 
extended component. The TC^C0)/TC^C0) ratio at 23.6'' 
resolution ranges from ~ 7 to ~ 1 1 when going from the in- 
ner to the outer part of the galaxy, in agreement with a mixture 
of dense and diffuse gas evolving from completely molecular 
in the inner galaxy to half atomic in the outer galaxy. 

Atomic thick disks are very common. Evidence for a 
molecular thick disk exists in at least the edge-on NGC 891 
galaxy. A thick molecular disk about 3 times as wide as the 
dense central CO layer was detected in the inner part of our 
Galaxy. We thus interpret the extended component of M5 1 
as a diffuse CO thick disk. The underlying physical inter- 
pretation of the ^^CO (1-0) emission is different. If the gas 
is dense, it fills a small fraction of the interstellar volume; it 
is confined by ram or turbulent pressure (if not gravitation- 
ally bound); and it is on the verge of forming stars. If the 
gas is diffuse, it is a warmer, low pressure medium filling a 
large fraction of the interstellar volume; it contributes more 
the mid-IR or PAH emission; and it is probably not gravition- 
ally bound or about to form stars. 
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APPENDIX 

A. CO LUMINOSITY AND MOLECULAR GAS 
MASSES 

A.l. Computations 

The CO luminosity, Leo, is estimated from the main beam 
temperature, r^b, as 



^^0 ^ 1 X] 7;ib(/,m,v) > A/mv, 



(Al) 



/,m,v 



where (/,m,v) are the indices of the position, position, veloc- 
ity data cube and A/mv is the volume of one pixel of this cube 



computed as 

A/mv = Av 
The uncertainty on this luminosity is computed as 



A/ Am/ D y 

\" \" V ^^7.6Mpcy 



SLco= l^r](l,m,v)STl^(l,m,v)Almv, 



(A2) 



(A3) 



where STj^\)(l^m,v) is the noise level for the pixel (l^m^v) and 
rjil^m^v) is the factor which reflects the correlations between 
pixels. If we assume that the noise level is independent of the 
velocity channel and that the velocity channels are uncorre- 
lated, we obtain 



^Lco = 



\ ny(l,m)r](l,m)ST^^(l,m)Almv, 



(A4) 



where n^il^m) is the number of velocity channels for the po- 
sition (l,m). The correlation between spatial pixels can be 
approximated as the inverse of the number of pixels in the 
resolution area, i.e.. 



r](l,m) = 



81og(2)A/Am 



(A5) 



where O and are respectively the FWHM major and minor 
axes of the Gaussian beam. It is independent of the pixel po- 
sition (/,m). 
The mass of the associated molecular gas, Mu , is then 



Mh ± SM^ = Xco (^co ± ^^co) , 



(A6) 



where Xco is the CO-to-H2 conversion factor. We will use the 
standard Galactic conversion factor 

Xco = 2.0 X 10^^H2cm"V(Kkms"^) (A7) 

= 4.35M0pc"V(Kkms"^). (A8) 

The last value includes a factor of 1.36 by mass for the pres- 
ence of helium. 

A. 2. Applications to M51 

The direct sum of the pixel brightness over the field 
of view observed with the IRAM- 30m map and between 
the [-200, +300] km s~^ velocity range indicates that the to- 
tal CO luminosity of M51 (including the companion) is 
1.442 X 10^ ±4 X 10^Kkms~^pc^. The same computation in 
the [-120, +120] km s~^ velocity range (where mainly M51a 
emits) gives 1.441 x 10^±3 x 10^Kkms"^pc^. Finally, the 
direct sum of the pixel brightness over the field of view where 
the line integrated emission is measured with a signal to noise 
ratio larger than 3 gives 1.428 x 10^ ±2 x 10^Kkms~^pc^. 
It is clear that the uncertainty on the result decreases when 
the summed volume of the cube is reduced because most of 
the original volume is devoid of signal. But reducing the vol- 
ume introduces some kind of bias. However, the total lumi- 
nosity does not vary significantly in both three results. We 
thus assume that most of the signal is contained within the 
smallest volume probed here. This in particular implies that 
the luminosity associated with the companion but outside the 
[-120,+120]kms~^ velocity range is negligible. The uncer- 
tainties are small compared to the absolute flux uncertainty 
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Table 7 

Total CO flux measured within the PdBI+30m data cube, as measured from 

integrated intensity images constructed using different techniques to isolate 

significant emission (see text). 



Method 



Parameters 



Total CO luminosity 
[xlQ^Kkms-^pc^l 



Direct sum 



[-120, +120] km s- 



9.1 



Sigma-clipping 


Tmb > Icr 


417 




Tmb > 2(7 


16 




Tmb > 3a 


7.0 




Tmh > 4cr 


4.6 




T^nb > 5a 


3.5 


Dilated mask 


a,p) = (3,1.2) 


10 




a,p) = (3,2) 


6.9 




a,p) = (4,1.2) 


7.6 




a,p) = (4,2) 


5.8 




(t,p) = (5,1.2) 


7.2 




a,p) = (5,2) 


5.6 


Smooth & mask 


(6',m) = (2.0,3) 


8.7 




(6',m) = (2.0,5) 


5.6 




(6',m) = (2.4,3) 


8.7 




(6',m) = (2.4,5) 


6.0 




(6',m) = (3.0,3) 


8.7 




(6',m) = (3.0,5) 


6.4 




(6',m) = (3.6,3) 


8.6 




(6>,m) = (3.6,5) 


6.5 



Hi prior 



lOkms"^ window 
20kms~^ window 
50 km s~^ window 



3.2 
5.0 
8.0 



Combined 



8.5 



(< 10%, Kramer et al. 20081 and the distance uncertainty 
(13%). We quote them once above to show that the radio- 
metric noise is negligible when estimating the CO luminos- 
ity. The surface where emission is detected is 38.3 square 
arcminutes, i.e., 1.9 x 10^ pc^. The associated total mass, 
mean brightness, and mass surface density are respectively 
6.2 X 10^ Mq, 7.6Kkms"\ and 33M0pc"2. These values are 
summarized in Table [T] 

The direct sum of the pixel brightness over the PAWS field 
of view and between the [-120, +120] km s~^ velocity range 
indicates that the final hybrid synthesis data cube contains 
a total CO luminosity of 9.1 x 10^Kkms~^pc^. This corre- 
sponds to a molecular gas mass of 4.0 x 10^ M©. The to- 
tal surface covered was 10.5 square arcminutes, i.e. 5.1 x 
lO^pc^. The mean brightness and mean surface density are 
respectively 1 8 Kkm s~^ and 77 M© pc~^ . 

B. MASKING TECHNIQUES 

In the previous section, we have seen that the line 0th or- 
der moment is much noisier when the velocity range used for 
integration includes a large number of channels devoid of sig- 
nal. This effect is amplified when computing the line 1st and 
2nd order moments. It is thus desirable to limit the velocity 
range of integration to channels where signal is detected. As 
the systematic motions are large for a galaxy, using a single 
velocity range for all the sky positions unavoidably implies 
moving the velocity range of integration from position to po- 
sition. Defining such moving velocity ranges is never straight- 
forward. In the case of the PAWS hybrid synthesis data set, it 
is even more complex because a large fraction of the flux lies 
at brightnesses barely detected (1 to 3cr). 

Moreover, masking out noisy channels always implies the 



risk to bias the result. We therefore explored several alterna- 
tive masking methods. In order to quantify the benefit over 
cost of different techniques, we compare the noise and signal 
aspects of 0th order moments images (Fig. 21 ) and the total 
luminosity found inside the kept velocity channels (Table [7]). 
Direct sum over a given velocity range produces an integrated 
intensity map that is dominated by noise (see Fig.[2TJa]) be- 
cause CO emission is typically only detected in a few velocity 
channels along each line-of- sight. But it produces a robust es- 
timate for the total CO luminosity. 

The first alternative we tried is a simple sigma-clipping 
method, whereby pixels containing emission with low signif- 
icance were excluded. We tested various brightness thresh- 
olds between 1 and 5 a, where a is the standard deviation 
of the noise fluctuations estimated using ~ 25 emission- 
free channels for each line-of- sight. The integrated inten- 
sity images were constructed by summing all unmasked (i.e. 
'good') pixels across the observed LSR velocity range, i.e. 
[-297.5, +297.5] km s~^ As expected, lower thresholds pro- 
duce maps that are dominated by noise peaks, especially in 
the interarm region. Higher thresholds produce maps with a 
cleaner appearance, but also do a poorer job of recovering the 
total luminosity (see Table [7]). For thresholds above ~ 3cr, 
it is evident that genuine low surface brightness emission in 
the interarm region is omitted. A map constructed using a 3cr 
threshold is shown as an example in Fig.[2TTb]. 

For the second method, which we call the dilated mask 
method, we defined islands of significant emission in the hy- 
brid synthesis data cube by selecting peaks above a thresh- 
old of t(j across two contiguous velocity channels. This pre- 
liminary mask was expanded to include all connected pix- 
els with emission greater than per. Several combinations 
of t and p values were tested, using t G [3,7] and p G 
[1.2,3]. The final mask was then applied to the original 
data cube, and the integrated intensity image was constructed 
by summing unmasked pixels within the LSR velocity range 
[-297.5, +297.5] km s~^ We consider the map obtained for 
(^, ^) = (5, 1 .2), shown in Fig.l2TTc], as the best map produced 
using this method. Higher t values tended to omit genuine low 
surface brightness emission, whereas lower t values include a 
larger number of spurious noise peaks, especially in the in- 
terarm region and at the edges of the survey field. Setting 
p=l.2 produced optimal results in the sense that faint emis- 
sion around the edge of bright structures is retained within the 
mask, but the algorithm does not misidentify large patches of 
noise as genuine emission. The total CO luminosity within 
the map shown in Fig. 
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Kkms ^pc^. 



c] is only 7.6 x 10^ 
however, indicating that a significant fraction of the CO emis- 
sion within the survey field is not recovered in this image. 

For the third method, which we call the smooth- and-mask 
method (e.g. [Heifer et al.[ 2003 ), we generated a mask by con- 
volving the original hybrid synthesis cube to a coarser spatial 
resolution using a Gaussian smoothing kernel with FWHM 6. 
The RMS noise for each sightline within the smoothed cube 
was estimated from ~ 25 emission-free channels, then pix- 
els in the smoothed cube with emission below a significance 
threshold mcr were blanked. After transferring this mask back 
to the original data cube, the integrated intensity image was 
constructed by summing unmasked pixels across the LSR 
velocity range [-297.5, +297.5] km s~^. As for the previous 
methods, we experimented with different combinations of the 
and m parameters. We consider the map obtained using a 
y.'6 kernel and a 5 a threshold to be the best produced using 
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Figure 21. Integrated intensity images constructed from the hybrid synthesis cube using different techniques: [a] direct sum of all pixels; [b] sigma-clipping 
method with a threshold of 5cr; [c] dilated mask method with (t,p) = (5, 1.2); [d] smooth-and-mask method with (0,m) = (3.6,5); [e] Hi prior method, with a 
50km s~^ integration window; [f] the combined method (see text). All images are presented on the same intensity scale, which is displayed near the bottom right 
image. 



this method. The total CO luminosity within this map, which 
is shown in Fig. 21 d], is 6.9 x 10^ Kkms~ 



-ipc^ 



For the fourth method, we integrated the PdBI+30m data 
cube over a narrow velocity range (which we refer to as the 
'integration window'), centered on the radial velocity at the 
peak of the Hi line profile for each line-of- sight. We call 
this method the 'Hi prior method'. The Hi velocity template 
was constructe d using the Hi data cube from THINGS ( Wal- 
ter et al.|2QQ8] ), which covers the entire disk of M51 at ~ 1 V^' 



resolution. Integration windows with velocity widths between 
10 and lOOkms"^ were tested. The rationale behind this ap- 
proach is that CO emission in nearby galaxies is mostly as- 
sociated with high brightness Hi emission (e.g.jSchrubaet al. 



[2011; Engargiola et al. 2003). Comparison with maps in the 
other panels of Fig. [21] suggests that CO emission in the nu- 
clear region of M5lTs not well-recovered by this method. 



One advantage of this approach, however, is that it yields an 
upper limit on the CO integrated intensity for pixels without 
detectable emission, and may therefore be more suitable for 
some types of quantitative analysis. The map for a 50 kms~^ 
integration window is shown in Fig.[2TTe]. The total luminos- 
ity in this map is 8.4 x 10^ Kkms~^ p?. 

Finally, we used the strengths of two of these methods to 
build our best mask. We started by building two 2D masks 
based on the second (dilated mask) and fourth (Hi prior) ap- 
proaches. In each case, we constructed a 3D dilated mask 
from the hybrid synthesis datacube. However, this technique 
also catches noise patches at "wrong" velocities. To remove 
these, we computed the associated centroid velocity map, and 
we then produced a 2D mask where the CO velocity centroids 
differed by less than 30kms~^ from a map of the Hi velocity 
field. Using this 2D mask, we could build a 3D filtered mask. 
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The difference between these two 3D filtered masks comes 
from the different thresholds used to produce initial 3D di- 
lated masks. The first dilated mask used (^,/?) = (4, 1), which 
allows us to find isolated faint point sources. However, these 
parameters also catches isolated noise at "forbidden" veloc- 
ities, which will considerably bias the line-of- sight CO cen- 
troid velocities. The comparison of the CO and HI centroid 
velocities may thus kill perfectly valid lines of sight. This 
happens in particular in the M51 central part. To avoid this 
effects, the second dilated mask used (t^p) = (10, 1.5). The 
high signal-to-noise ratio (10) starting point ensures that no 
outlier velocities will be caught in the 3D mask but it will 
miss isolated low brightness signal. The two 3d filtered masks 
thus identify complementary kind of signal. We finally make 
the union of these 3D filtered masks to obtain the final mask. 
The resulting 0-moment map, which has a total luminosity of 
8.5 X 10^Kkms~^pc^, is shown in Fig. 21 f]. This map con- 



tains 93% of the total luminosity computed from direct sum- 
mation. The comparison of these maps with the ones of the 
other methods shows that it is the most successful at recover- 
ing all genuine emission within the data cube. This mask was 
thus used to compute the 1st and 2nd order moments of the 
hybrid synthesis cube at V^ resolution. 

C MODELING THE IMPACT OF THE IRAM-30M 
ERROR BEAM 

In this appendix, we assess how much the error beams of 
the IRAM-30m contribute to the extended emission. To do 
this, we first convolved a model of the source emission in the 
PAWS field of view 1) with the ideal 30m beam and 2) with 
a model of the true 30m beam, including the diffraction pat- 
tern and the error beams. We then subtracted the convolved 
maps to estimate the contribution of the diffraction pattern 
and error beams to the extended emission. We first describe 
the 30m beam and the source models used, before discussing 
the results. 

C.l. Beam model 

After measuring the 30m beam pattern at 3.4, 2, 1.3 and 
0.8 mm from total power sc ans across th e Moon around full 
Moon and new Moon, Grev e et al.| ( |1998| ) modelled the beam 
using antenna tolerance theory. They deduced that the beam 
consists of the diffracted beam, two underlying error beams, 
which can be explained by the panel dimensions, and a beam 
deformation mostly due to large-scale transient residual defor- 
mations of the telescope structure. We scaled their results to 
the frequency of the ^^CO (1-0) line to model the error beams 
we use. The diffraction pattern was computed for an antenna 
illuminated by a Gaussian beam of 12.5 dB edge taper and a 
ratio of the secondary-to-primary diameter (blockage factor) 
of 0.067, using the prescription of Goldsmith (1998). Fig- 
ured 



; [22] shows the properties of the resulting beam and Table [8] 
listsme typical scales, amplitudes, and contributions to the 
total power of the different beam components. 

The diffraction pattern and two of the three error beams are 
linked to the structure of the primary mirror of the 30m, which 
have not changed since the study of Greve etal.| ( |199 8). Hence 
the typical angular scales of the different error beams stay 
constant over time. However, the 30m primary surface ac- 
curacy improved in 2000 after several holography campaigns. 
Moreover, the thermal balance of the telescope primary sur- 
face was also improved around this time. Both effects ex- 
plain why the beam efficiencies of the 30m are significantly 
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Figure 22. Characteristics of the IRAM-30m beam at the frequency of the 
^^CO (1-0) line as a function of the opening angle in linear (left column) and 
logarithmic (right column) scales. Top: Cut of the beam profile, normalized 
to unity at ^ = 0^^ The black curve is a combination of the diffraction pattern 
and the three Gaussian error beams. The green, blue, and red curves displays 
the three error beams. Middle: Beam power contained in circular annuli, 
normalized so that the integral is unity. Bottom: Beam power integrated over 
the solid angle sustained by the opening angle, normalized to 100%. 



higher today than the values provided by |Greve et al.| ( p^98] ): 
e.g., 0.78 instead of 0.72 at the ^^CO (1-0) frequency. We 
used slightly different values of the 30m beam efficiency in 
our analysis for historical reason. As no newer detailed mea- 
surements of the error beams are available, we used 0.72 for 
consistency with the work by Greve etal.| ( |1998| ) to model the 
beam here, while we applied 0.75 to the data. As the cur- 
rent best estimate of the beam efficiency is 0.78, both values 
present only variations of a few percent and do not alter our 
conclusions. In particular, our estimates of the error beam 
contribution can be seen as a firm upper limit. 

C.2. Source model 

The best proxy for the source model is the result of this 
project, i.e. the distribution of the emission measured in the 
PAWS field of view. We made two different source models. 
First, we want to check whether the error beam contribution 
of the compact sources alone can account for the extended 
emission. We started with the \" PdBI-only data cube, and we 
set to zero all pixels whose brightness was below 3 times the 
noise level in order to avoid spurious sources. We multiplied 
all pixel brightnesses by a factor of about 2, needed to recover 
the total flux in the PAWS field of view. This assumes that 
all extended emission is spurious and that any flux associated 
with it should indeed belong to compact sources. We will 
refer to this cube as the compact model. 

In contrast, we assume for the second model that the ex- 
tended emission is genuine. To model the extended compo- 
nent, we start with the 6" subtracted cube, in which we set to 
zero all pixels whose brightness was below 3 times the noise 
level. We then added all the brightness of the \" PdBI-only 
data cube whose brightness is above 3 times the noise level. 
The flux contained in this model is then 3% larger than the 
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Table 8 



Parameters used to model the IRAM-30m beam at 115.271202 GHz. Adapted from |Greve et al.| ( |1998) 
Main beam 1st error beam 2nd error beam 3rd error beam 



Origin 

Correlation length 
Surface rms 
Beamwidth (FWHM) 
Power amplitude 
Relative power 



Diffraction 
pattern 



21.3'' 

1.00 

71.8% 



Large scale 

deformations 

2.5-3.5m 

< 0.06mm 

230'' 
9.0 X 10-4 

7.4% 



Panel Frame 

misalignments 

1.5-2.0m 

0.055 mm 

330" 

2.7 X 10-4 

4.7% 



Panel 

deformations 

0.3-0.5m 

0.055 mm 

1840" 
3.1 X 10-5 

16.1% 
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Figure 23. Comparison of the spatial distribution (from top to bottom) of the peak intensity, integrated intensity, centroid velocity, and the line full width at 
half maximum {i.e., 2.35 times the standard deviation in velocity) of the ^^CO (1-0) emission from the actual observations (left column), and the two models 
convolved with the modeled 30m beam (middle and right column). The angular resolution is indicated by a circle in the bottom left corner of each panel. The 
intensity scale is shown on the right-hand side of each panel. The 3 images of each row share the same intensity scale to facilitate a meaningful visual comparison. 
The crosses on the images show the positions of the spectra displayed in Fig.|24] Other plot annotations are the same as in Fig. [13] 



flux measured in the PAWS field of view. We will refer to this 
cube as the summed model. 

C.3. Modeling results 

We convolved each model with the true 30m beam model. 
We normalized the result by the main beam efficiency to get 
to the main beam temperature scale. Fig. [23] compares the 
moments of both modeled 30m emission witn the measured 
emission. All the moments were computed using the same 
position-position- velocity mask (the one of the data) to en- 
sure a meaningful comparison. The compact model produces 
much higher peak temperatures and integrated emission in the 
arms, while it understimates the flux in the inter-arm region. 



The summed model provides a much better match to the ob- 
servations. 

We subtracted from the previous cubes the associated 
brightness model convolved with the ideal 30m beam. This 
yields the contribution of the error beams to the signal mea- 
sured with the 30m. In particular, it allows us to estimate 
the error beam contribution to the extended emission, as we 
showed in Section [3X4] that the flux of the extended emission 
is mostly structured at scales larger than 36'^ i.e. at scales 
larger than the ideal 30m beam FWHM. For both source mod- 
els, the flux scattered into the error beams is less than 20% the 
total flux in the PAWS field of view but the peak brightness 
due to the error beam contribution is only 55 mK, i.e., about 
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Figure 24. Comparison of the spectra of our summed model (in black), the extended emission as measured in the 6" subtracted cube (in red), and the error beam 
contribution (in green), at six positions where the extended emission shows different characteristics. These positions are displayed as crosses on the images of 
the extended emission moments in Fig.|13| The bottom row shows a brightness zoom of the top row. 
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Figure 25. Comparison of the spatial distribution of the line full width at half 
maximum (i.e., 2.35 times the standard deviation in velocity) of the ^^CO 
(1-0) extended emission measured at 6" . Different masking techniques were 
used to test the impact of the error beams on the computation of the line 
width (see text). The angular resolution is indicated by a circle in the bottom 
left corner of each panel. The intensity scale is shown on the right-hand 
side of each panel. The images share the same intensity scale to facilitate 
comparison. The major and minor axes are displayed as perpendicular dotted 
lines. The dotted circles show the two inner corotation resonances at radii 
equal to 23^^ and 55^^ while the d ashed circle sh ows the start of the material 
arms at a radius equal to 85^^ ^Meidt et al.|2013) . 

3.5 times the median noise level of the 30m observations. The 
typical angular scales of the error beams are large, implying 
that the flux is scattered at low brightness level over wide re- 
gions of the sky. The median value of the error beam contri- 
bution is 19mK, while pixels brighter than 5 times the noise 
level have a median of 140 mK in the extended emission mea- 



sured at &' (see Sect. |3.2.T] ). Moreover, flux above O.IK in 
the extended emission cube represents 79% of the flux in this 
cube (we subtracted the flux contributed by the error beams in 
the same mask to compute this value). We deduce that the er- 
ror beams contribute less than 20% of the flux in the extended 
emission, i.e., less than 10% of the total flux in the PAWS field 
of view. We speculate that the baselining of the 30m spectra 
removed a large fraction of the flux associated with the error 
beams. 



Fig. 24 compares the spectra of our summed model, the ex- 
tendecTemission, and the contribution from the error beams 
at six positions, which sample different ratios of compact to 
extended brightness as well as different characteristics of the 
extended emission. These spectra illustrate that the emission 



associated with the error beams has a very different signature 
in space and velocity from the extended emission distribu- 
tion. In particular, the spectra resulting from the error beams 
have much wider linewidths than the observed spectra of the 
extended emission. We finally investigate the impact of the 
error beams on the moments of the extended emission. We 
recomputed the moments of the 6^' extended emission cube in 
two different ways. First, we only include the pixels brighter 
than two times the peak brightness (/. ^. , 0. 1 K) due to the error 
beams. In the second test, we include only the pixels brighter 
than four times the contribution of the error beam for each 
source model. The three resulting moment maps are similar 
to the original one. This is also true for the second moment 
or linewidth displayed in Fig. [25] We thus conclude that the 
large linewidths of the extendeoemission are genuine. 

D. ESTIMATING THE VERTICAL VELOCITY 
DISPERSION 

The velocity of a parcel of gas in the galaxy frame is 

Vgal = (Vp,Vz), (Dl) 

where (Vp,Vz) respectively are the in-plane and perpendicular 
to the galactic plane velocity components. In the observing 
frame, the Doppler effect allows us to infer an observed ve- 
locity, Vobs, which is the projection of Vgai along the line of 
sight. If Up is the projection of Vp on the galaxy major axis 
and 1 is the inclination of the galaxy plane onto the line of 
sight, we get 



Vobs = i^p sinB + v^ cosB. 



(D2) 



In the optically thin limit, a line is an histogram of all the 
Vobs values along each line of sight. For a resolved perfectly 
face-on galaxy, 

Vobs = v^. (D3) 

The centroid velocity would thus be constant and the 
linewidth would give the dispersion of the velocity distribu- 
tion along the vertical axis under the condition that the veloc- 
ity has a symmetric distribution of velocity along any line of 
sight. For a resolved edge-on galaxy. 



^obs ~ ^p- 



(D4) 



If the velocity is only rotational and the density is much higher 
in the spiral arm, the centroid velocity is then biased towards 
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the velocities in the spiral arm and the velocity dispersion 
gives an idea of the velocity content along the line of sight. 
This is why we still have some information about the rotation 
curve of the galaxy. 
For any other inclination, the centroid velocity, Vcent. is 



Vcent = (Vobs) = {Up) sinB, 



(D5) 



as long as the galaxy disk is not warped and the distribution of 
the velocity perpendicular to the galactic plane is symmetric, 
i.e. {vz) = 0. Now, the velocity dispersion is computed as 



<(v, 



)2) = / [v^ cosB+ {up-{up)) sinB]^) . (D6) 



obs " 



Assuming that there are no correlations between vertical and 
horizontal motion, we obtain 

((Vobs -Vcent)') = {v]) cos^ B+ ^ (w^ - (w^))') sin^B. (D7) 

If we also assume that the velocity field is only made of a 
systematic motion (Wsys) parallel to the galaxy plane plus an 
isotropic 3D turbulent motion of typical dispersion (aturb), 
then 

{vI) = cjI,^ and (D8) 



(Wp-(Wp))'^=Crt^urb+((^sys-(Wsys))^). 



(D9) 



Note that (wsys) = (^/?) ^iid the turbulent component yield only 
one time dturb because Up is an unidimensional velocity com- 
ponent. Hence, 

((Vobs-Vcent)') =Cr^2^rb + ([(^sys-(Wsys)) sinBj'), (DIO) 



where / [(wsys - (wsys) sinB)] \ is the contribution of the sys- 
tematic in-plane motions inside the beam to the line 2nd order 
moment. We can estimate this contribution as 



[(Wsys-(Wsys)) Sln B] ^ 



|grad(Vcent)| 



V^Sbl. 



(DID 



where 6 is the resolution beamwidth of the observations 
and |grad(Vcent)| is the modulus of the centroid velocity 2D- 
gradient. In our case, these last two quantities can be mea- 
sured. This thus yields an estimation of the galaxy vertical 
velocity dispersion as the square root of 



(^z> = ^turb ^ ((Vobs -Vcent)' )- 



n2 



|grad(Vc, 



2.35 



. (D12) 



E. ADDITIONAL MATERIAL 



This section displays additional material that completes the 
description of the observations and the presentation of the data 
cubes. All these tables and figures are only made available 
online. 

Table [9] presents the weather conditions during the IRAM- 
30m observing run. Table 10 exhaustively lists the interfer- 
ometric sessions observed auring the IRAM large program 
PAWS. Table 1 1 gives the evolution of the calibrator fluxes as 
a function of time, a key intermediate in the absolute ampli- 
tude calibration of the data. 



Eigures [26| and [27] present the channel maps of the IRAM- 
30m single-dish data for the ^^CO and ^^CO (1-0) lines. Eig- 
ures 28 and 28 show the channel maps of the hybrid synthe- 
sis data cube at \" angular resolution, the main data product 



Table 9 

Detailed parameters for the IRAM-30m observations 



Observing date 


Time^ 


T"- b 

^sys 


Water vapor 




hours 


K[r;] 


mm 


18/05/2010 


3.8/8.0 


279 


4.3-8.8 


19/05/2010 


2.9/8.0 


275 


3.6-6.1 


20/05/2010 


3.4/8.0 


273 


2.6-8.5 


21/05/2010 


3.6/9.0 


287 


5.8-9.0 


22/05/2010 


3.6/8.0 


297 


6.5-9.4 



^ Two values are given for the integration time: the on-source time and the 

telescope time. 

^ The Tsys values are given for the ^^CO (1-0) frequency. 



Table 10 

Detailed parameters for the PdBI observations. 



Config. A^ar 



Mosaic 
# 



Int. Time^ 
hours 



Seeing 
arcsec 



Obs. date 



D 


5 


Top 


2.0/09.5 


150-400 


1.70 


29 Aug. 2009 




5 


Bot 


2.2/09.5 


200-400 


1.70 


02 Oct. 2009 


C 


6 


Top 


4.3/09.5 


180-250 


0.68 


27 Oct. 2009 




6 


Bot 


3.3/08.5 


180-300 


0.85 


28 Oct. 2009 


B 


6 


Top 


1.8/08.0 


200-250 


0.37 


20 Feb. 2010 




6 


Top 


1.2/05.0 


200-250 


0.60 


02 Mar. 2010 




6 


Top 


1.5/05.0 


160-230 


0.50 


04 Mar. 2010 




6 


Top 


3.3/07.5 


180-280 


0.29 


11 Mar. 2010 




6 


Bot 


5.1/11.0 


150-250 


0.34 


08 Mar. 2010 




6 


Bot 


4.2/11.5 


200-300 


0.50 


09 Mar. 2010 




6 


Bot 


4.6/09.0 


160-220 


0.30 


12 Mar. 2010 


A 


6 


Top 


1.3/03.0 


180-200 


0.32 


14 Dec. 2009 




6 


Top 


3.0/08.5 


180-220 


0.28 


17 Dec. 2009 




6 


Top 


5.2/11.0 


160-230 


0.15 


18 Jan. 2010 




6 


Top 


2.2/05.0 


200-600 


0.27 


22 Jan. 2010 




6 


Top 


3.6/07.5 


180-230 


0.27 


12 Feb. 2010 




6 


Top 


1.5/04.5 


160-230 


0.16 


13 Feb. 2010 




6 


Bot 


3.2/05.5 


170-220 


0.31 


15 Dec. 2009 




6 


Bot 


1.5/07.5 


180-400 


0.35 


03 Jan. 2010 




6 


Bot 


3.0/08.5 


180-400 


0.44 


05 Jan. 2010 




6 


Bot 


3.2/05.5 


180-230 


0.18 


23 Jan. 2010 




6 


Bot 


2.3/06.5 


250-300 


0.41 


29 Jan. 2010 



^ Two values are given for the integration time: the good on-source time 
if observed with 6 antennas) and the telescope time. 



(as 



of the PAWS large program. Figures [30| and [31] overlay the 
channel maps of the bright compact emission as contours over 
the channel maps of the faint extended emission displayed in 
color. 



Finally, Figures 32 33 34 and 35 summarize the properties 
of the decomposition orthenybriasynthesis emission into the 
PdBI-only and the subtracted emissions, respectively at 1, 3, 
6,and22.5'^ 



PAWS REVEALS A THICK DISK OF DIFFUSE MOLECULAR GAS IN THE M51 GALAXY 31 

Table 11 

Flux (in Jy) of the amplitude and phase calibrators used during the PdBI 
calibration. 



Date 


1418+546 


1308+326 


Jl 332+473 


29-aug-2009 


1.05 


— 


0.46 


02-oct-2009 


1.10 


2.58 


— 


27-oct-2009 


1.06 


— 


0.22 


28-oct-2009 


1.02 


2.58 


— 


14-dec-2009 


1.11 


— 


0.22 


15-dec-2009 


1.10 


2.02 


— 


17-dec-2009 


1.16 


— 


0.21 


03-jan-2010 


1.23 


2.26 


— 


05-jan-2010 


1.23 


2.16 


— 


18-jan-2010 


1.22 


— 


0.23 


22-jan-2010 


1.16 


— 


0.24 


23-jan-2010 


1.16 


1.93 


— 


25-jan-2010 


1.22 


1.83 


— 


29-jan-2010 


1.06 


1.83 


— 


12-feb-2010 


0.98 


1.91 


0.22 


13-feb-2010 


1.00 


— 


0.22 


20-feb-2010 


0.92 


— 


0.18 


02-mar-2010 


0.92 


— 


0.19 


04-mar-2010 


0.94 


— 


0.20 


08-mar-2010 


0.94 


1.92 


— 


09-mar-2010 


0.96 


1.95 


— 


ll-mar-2010 


0.96 


— 


— 


12-mar-2010 


0.96 


— 


— 
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Figure 26. Channel maps of the ^^CO (1-0) emission obtained with the IRAM-30m telescope. The velocity in kms~^ of each channel is displayed in the top left 
corner of each panel. The intensity scale (in r^b) is shared by all the panels and it is displayed in the bottom right corner of the figure. The major and minor axes 
are displayed as perpendicular dotted lines. The dotted circles show the two inner corotation resonances at radii equal to 23^^ and 55^^ while the dashed circle 
shows the start of the material arms at a radius equal to 85^^ (Meidt et al. 2013). 
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Figure 27. Channel maps of the ^^CO (1-0) emission obtained with the IRAM-30m telescope. The figure layout is the same as for Fig. [261 
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Figure 28. Channel maps of the ^^CO (1-0) emission obtained from the combination of IRAM-30m and IRAM-PdBI observations. The velocity of each channel 
in kms~^ is displayed in the top left corner of each panel. Only the negative velocity channels are shown here. The intensity scale (in r^b) is shared by all the 
panels and it is displayed in the bottom right corner of the figure. This intensity scale is saturated to emphasize the pixels with significant emission. The major 
and minor axes are displayed as perpendicular dotted lines. The dotted cir cles show the two inner corotation resonances at radii equal to 23^^ and 55^^ while the 
dashed circle shows the start of the material arms at a radius equal to 85 '^ (Meidt et al.|2013| . 
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Figure 29. Same as Fig.l28]but for the positive velocity channels. 
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Figure 30. Contours of the channel maps of the signal-to-noise ratio of the co mbin ed PdBI-i-30m cube laid over the channel maps of the resolved emission. The 
contour levels start at a signel-to-noise ratio of 8. The figure layout is as in Fig. [28] 
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Figure 31. Same as Fig.jSOlbut for the positive velocity channels. 
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Figure 32. Comparison of the spatial distribution at ~ 1^^ -resolution (from top to bottom) of the peak intensity, integrated intensity, centroid velocity, the line 
full width at half maximum (i.e., 2.35 times the standard deviation in velocity) and, rms noise of the ^^CO (1-0) emission for the hybrid synthesis (PdBI + 30m, 
left column), the PdBI-only (middle column) and the subtraction of the PdBI-only from the hybrid synthesis cubes (right column). The angular resolution is 
shown as a circle in the bottom left corner of each panel. The intensity scale is shown on the right of each panel. The major and minor axes are displayed as 
perpendicular dotted lines. The dotted circles show the two inner corotation resonances at radii equal to 23^^ and 55^^ while the dashed circle shows the start of 
the material arms at a radius equal to 85^^ ^Meidt et al.|2013| . 
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Figure 33. Same as Fig.[32]but the comparison is done here at an angular resolution of 3" 
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Figure 34. Same as Fig.[32]but the comparison is done here at an angular resolution of 6" 
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Figure 35. Comparison of the spatial distributions at 22. 5^^ -resolution of the peak intensity (top), integrated intensity, centroid velocity, 2.35 times the standard 
deviation and rms noise (bottom) of the ^^CO (1-0) emission for IRAM-30m cube (left column), the (i" PdBI-only data cube smoothed at 22.5^^ (right column), 
and the subtraction of the two previous cubes (middle column). The angular resolution is displayed as a circle in the bottom left corner of each panel. The 
intensity scale is shown on the right side of each panel. The major and minor axes are displayed as perpendicular dotted lines. The dotted circles show the two 
inner corotation resonances at radii equal to 23^^ and 55^^ while the dashed circle shows the start of the material arms at a radius equal to 85^^ ( Meidt et al.,2()13) . 



